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CHAPTER I 


INTRODUCTION 

The purpose of this study is to develop a methodology to estimate the reliabil- 
ity of engineering systems that encompass several disciplines. The methodology is 
implemented using the NESSUS probabilistic analysis code, which has mostly been 
applied exclusively in the discipline of structural engineering. In order to apply 
the NESSUS probablistic structural analysis code to analyze a multi-disciplinary 
engineering system, the equivalences between system behavior models in different 
disciplines are investigated, and the effect of physical interaction among the failure 
modes is quantified in this study. 

System reliability analysis is a method of estimating the effects of uncertainties 
in an engineering system on the probability of successful performance. Usually, 
an engineering system consists of multiple subsystems and components, which may 
require the knowledge of different disciplines of engineering. Such disciplines may 
include structural engineering, mechanical engineering, heat transfer theory, fluid 
mechanics, electrical engineering, etc. Such a system is called a multi-disciplinary 
engineering system. The reliability analysis of any engineering system usually begins 
with the identification and reliability computation of individual failure modes within 
the system. Then the reliability analysis of the overall system can be carried out. 

Traditionally, reliability methods have primarily concentrated on failures in one 
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particular discipline, e.g. structural analysis, not on an overall system which con- 
sists of multiple disciplines. Furthermore, conventional methods of system reliability 
estimation usually only consider the statistical correlation between individual failure 
events, ignoring the fact that more often than not, those individual failure modes also 
have a physical correlation. This leads to inaccuracy in system reliability estimation. 

The method presented in this report uniquely computes the failure probability 
interactions between different modes and overal system failure probability through 
the imposition of one failure mode on another field and reanalysis of the latter. This 
method is used to compute the probabilities of critical system failure events after 
accounting for the contributing non-critical failure modes in all different fields. How- 
ever, it is not an easy task to estimate the reliability interactions between different 
failure modes. The success of such a method primarily depends on the availability of 
effective reliability tools. The software system NESSUS developed under the leader- 
ship of NASA Lewis Research Center is uniquely suited for this purpose. Currently, 
this code has been applied primarily to the structural engineering problems. In or- 
der to perform system reliability analysis including the interactive failure modes, this 
study uniquely develops behavior analogies between the structural model and heat 
transfer model, and between the structural model and fluid mechanics model. By 
doing so, the probability estimation of heat transfer and fluid mechanics failures can 
be pursued similarly to structural reliability analysis. 

The objective of this research project is to develop a method, using system relia- 
bility theory, for the reliability estimation of multi-disciplinary engineering systems. 
The method is implemented on the software system NESSUS (Numerical Evalu- 
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ation of Stochastic Structures under Stress) developed by NASA Lewis Research 
Center. An example application to a three-discipline system involving mechanical 
stress-strain behavior, heat transfer and fluid mechanics is provided. In order to 
compute the individual failure mode probability of non-structural problems such as 
heat transfer and fluid mechanics within NESSUS, it is necessary to develop a new 
methodology for the analysis of heat transfer problems using the concept of equiva- 
lence between the computational quantities in structural analysis, such as stiffness, 
displacement vector, load vector, etc. and similar quantities such as conduction, 
temperature distribution and heat flux in heat transfer theory, flow velocity, pres- 
sure and flow factor in fluid mechanics. This is the first important contribution of 
this study. ' 

The second important contribution is the method for the computation of the 
physical dependence of critical failure mode probabilities on non-critical failure modes 
in various disciplines. This involves the imposition of the non-critical modes and 
reanalysis of the system with appropriate discipline equivalences, for various levels of 
progressive damage. The combination of these two ideas - inter-disciplinary analogies 
and physical failure mode correlation - makes a reliability analysis program such as 
NESSUS very powerful for application to a variety of multi-disciplinary systems. 

The concepts and methods discussed above are examined in detail in the next 
four chapters of this report. In Chapter II, the basic reliability analysis concepts 
for individual component-level and system-level events are reviewed, and their im- 
plementation in the NESSUS program is described. In Chapter III, the behavior 
analogies between the structural analysis model and heat transfer problem, and be- 
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tween the structural analysis model and fluid mechanics model are developed. The 
finite element numerical examples with NESSUS/FEM are demonstrated for this con- 
cept. Chapter IV consists of two major parts: in the first part, the failure probability 
analyses for individual events including structural, heat transfer and fluid flow failure 
modes are performed using NESSUS; in the second part, the system failure proba- 
bility is studied. The effect of non-critical failure events of heat transfer and fluid 
mechanics upon a critical structural failure event is investigated, followed by system 
reliability analysis with the consideration of physically correlated component-level 
events. A numerical example of system reliability analysis of a multi-disciplinary 
system consisting of structural, heat transfer and fluid mechanical modes is demon- 
strated. The conclusions and recommendations of the study are summarized in 
Chapter V. 
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CHAPTER II 


SYSTEM RELIABILITY ANALYSIS 

Individual failure modes and effects 

An engineering system consists of a number of functional components. Before the 
system- level analysis begins, the modes of failure for individual components should 
be specified. The analyses of the failure modes and effects can be carried out by 
starting at the component level and expanding upward to the whole system. A 
failure mode is the manner by which a failure is observed. All units in a system are 
designed to fulfill one or more functions. A failure is thus defined as non-fulfillment of 
one of these functions. Analytically, each failure mode has a corresponding limit sate 
which separates the design space into “failure” and “safe” regions. The probability 
of failure, Pf, is denoted as 

Pf = P[g < 0] (i) 

where g is the value of the performance function g(X). The limit-state is denoted 
by the equation g(X) — 0. 

An exact solution of Pf can be obtained by the integration of the multiple integral 
denoted as 

Pf = l fx(x)dx (-) 

~:{X)S0 
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where f(X) is the joint probability density function of the vector of uncertain vari- 
ables X. 

In general, the solution of this multiple integral is too complicated to obtain. 
This is not only because the individual distributions are not always available but 
also because the integral is multi-dimensional for a realistic problem and is difficult 
to evaluate. Therefore, for practical purposes, efficient approximate analysis tools 
are needed. 

Fig. 1 illustrates the concept of the first-order approximation to the limit state 
for an estimate of the failure probability. 


The uncertain variables [X ) are all transformed to equivalent un correlated stan- 
dard normal variables (u). The most probable point MPP of the limit state is defined 
at the minimum distance /? from the origin to the limit state surface. Therefore, the 
first-order estimate of the failure probability is 

P } = $(-/?) 

where $ is the distribution function of a standard normal variables. 

In the NESSUS computer code, this is referred to as the Fast Probability Inte- 
gration (F PI) method. The limit state is constructed as; 

g=Z(X)-Z o = 0 m 
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Figure 1: Failure Probability Estimation 
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where Z 0 is a real value of the random variable Z(X ), which is a performance function 
or a response function, such as stress, displacement, temperature, etc. 

Z(X) = Z(Xi,X 2 , ..., X n ) (5) 

where X{(i = 1,2, ...,n) are the input random variables. 

The NESSUS program searches for the MPP by computing the sensitivities of the 
limit state to the random variables using iterative perturbation (in NESSUS/PFEM). 
and using these sensitivities to obtain a mean value first order (MV FO) or second- 
order (MV SO) estimate of the failure probability (in NESSUS/FPI). 

By using a first or second-order Taylor’s series expansion around the MPP, u*, 
the exact <7(u)-function is replaced by the first-order polynomial, gi(u), 

71 

gi(u) = a 0 -I- ^2 a,i(ui — u *) /g\ 

i=i v 

ox a second-order polynomial, g 2 (u), 

g 2 (u) = a 0 + J2 a x (m -u*) + ]T &<(u f - - <) 2 + ^ £ <?,•>(«,• - u*)(u 3 - u]) (7) 

‘- 1 «=i *=i i=i 

where the coefficients can be obtained by perturbation. 

Once these functions are obtained, the MPP is found. The probability of failure 
can be computed easily using Equation 3. This is the mean value first order (MV FO) 
estimate of the failure probability. This is improved using the Advanced Mean Value 
(AM \ ) analysis. Point probability estimate is made using specific limits for Zo, and 
the cumulative distribution function (CDF) is obtained by varying Zq. 
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Figure 2: Series System Structure 

NESSUS/FEM employs innovative finite element technology and solution strate- 
gies. It provides a choice of algorithms for the solution of static and dynamic prob- 
lems, both linear and nonlinear, together with an interactive perturbation analysis 
algorithm to evaluate the sensitivity of the response to small variations in one or 
more user-defined random parameters. 

NESSUS/FPI (Fast Probability Integrator) is used to evaluate structural response 
cumulative distribution functions (CDF). There are two methods in the code, the 
first-order reliability method and the advanced first-order reliability method. In gen- 
eral, the structural performance or response functions (e.g., stresses, displacements, 
vibration frequencies) are implicitly defined and each function evaluation may require 
intensive computation. The AMVFO ( Advanced Mean Value First Order) method 
reduces the computational burden and is the main probabilistic tool in NESSUS. 
NESSUS/PFEM automates the AMV procedure by integrating the FPI code and 
the FEM code. 
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System and component-level failure modes 

After the individual reliability analysis is completed, one then proceeds to the 
system or subsystem level analysis. System failure may occur due to a combination 
of any of the individual component failure modes. Many physical systems that are 
composed of multiple components can be classified as series-connected or parallel- 
connected systems, or combinations of series and parallel conditions. Description of 
these simple system structures is as follows. 

• Series System 

A system that is functioning if and only if all of its n components are functioning 
is called a series system structure. Fig. 2 illustrates such a system. 

If E{ denotes the failure mode i, then the failure of a series system is the event 

Ef = E x U E 2 U ... U E n (8) 

Then the failure probability of the system is 
Pi = P(Ej) (9) 


If each failure mode Ei is represented by a limit state = 0 in basic variable 

space, the failure probability can be obtained by the integration denoted as 


Pi= L x -I h(x)d * 


( 10 ) 
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s region : g(x) < 0 


in two dimensions 


a 



b 


Figure 4: Parallel System Structure 

where X represents the vector of all the basic random variable (loads, material 
properties, etc.) and is the domain in X defining failure of the system. 

This is defined in terms of the various failure modes as 0 t (AT) < 0. In two- 
dimensional X space, expression (10) is defined in Fig. 3. 

• Parallel System 

A system that is functioning if at least one of its n components is functioning is 
called a parallel system structure. A parallel structure of order n is illustrated in 
Fig. 4. 
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Figure 5: Combined System Structure 
In this case, the system failure event can be written 

Ej = Ei n e 2 n ... n E n (n) 

• Combined Series-Parallel System Structure 

This refers to systems which are a combination of series and parallel structures. 
Fig. 5 shows an example of such systems. 

The failure event of this system is written, for example, as 

E f = [Ei n (E 2 U E 3 )} U E 4 (12) 

It should be noted that not all engineering systems can be represented simply as 
described above. Practical systems may be more complex and need more effort to 
model. 
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System reliability computation 

In NESSUS, two methods have been implemented for system reliability compu- 
tation [4]: (1) probabilistic fault tree analysis combined with importance sampling 
(Torng et al, 1992) and (2) a structural reanalysis procedure to accurately estimate 
the failure regions for various critical failure modes affected by progressive damage 
(Mahadevan et al, 1992). 

Consider an engineering system subject to a sequence of loads (duty cycles) and 
which may fail in any one (or more) of a number of possible failure modes under any 
one load in the loading sequence. The total probability of the system failure may 
then be expressed in terms of the individual mode failure probabilities as 

p j = P{E X ) u p(e 2 n s x ) u p(e 3 n s 2 n s x ) u P(e 4 n s 3 n s 2 n Si) u ... (13) 

where E t denotes the “failure of the system due to failure in ith mode and Si denotes 
the complementary survival event of the ith mode. 

Since P(E 2 D Si) = P(E 2 ) - P(E 2 n E x ), Eq. 13 may be written also as 

Pj = P(E 1 ) + P(E 2 ) - P(E X n e 2 ) + P(E 3 ) - p(e x n e 3 ) 

—P(E 2 n e 3 ) + P(E X n e 2 n e 3 ) -j- ... 

where (E x D E 2 ) is the event that failure occurs in both modes 1 and 2, etc. 

Since it is not always an easy task to determine the joint probabilities of more 

than two failure modes, the following approximation methods can be used to predict 
the system reliabilities. 
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• First-order bounds 


The probability of failure for the system can be expressed as Pf = 1 — P{S), 
where P(S) is the probability of survival. For independent failure modes, P(S) can 
be represented by the product of the mode survival probabilities, or, noting that 
P(Si) = 1 - P(E t ), by 

p/ = i - ni 1 - < 15 > 

»=i 

where, as before, Pf is the probability of failure in mode i. This result can be shown 
to be identical with Eq. 14. It follows directly from Eq. 14 that, if P(E{) <!C 1, then 
Eq.15 can be approximated by [Freudenthal et al., 1966] 

P/ = E P(Bi) < 16 > 

i=l 

In the case where all failure modes are fully dependent, it follows directly that the 
weakest failure mode will always govern system failure, irrespective of the random 
nature of the strength. Hence 

Pf = max[P(FJ,)] (17) 

1 = 1 

Equations 15 or 16 and 17 can be used to define relatively crude bounds on 
the failure probability of any system of the series types when the failure modes are 
neither completely independent nor fully dependent. These are Cornell’s first-order 
bounds: 

max P(Ei) < P(U? =l Ei) < £ P{Ei) (18) 
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• Second-order bounds 


For some practical systems, the above first-order bounds may be too wide to 
be meaningful. For more accurate estimation, second-order bounds have been de- 
veloped. There are various second-order bounds in the literature [Kounias, 1968; 
Vanmarcke, 1973; Hunter, 1976; Ditlevsen, 1979]. Cruse et al (1992) derived second- 
order bounds which are independent of any ranking of the failure events[l]. The 
upper bound is 

p ( u -i ^ dt P ( E i) - max E m * x P(EiEj), max E P(EiEj)]} (19) 

'=2 J< ‘ --i=l 

The lower bound is 


n max(i,j) 

P(ur =1 £ t ) > m^X n {P(Ej) + max[(P(^)- E P{E { E k )); 0]} (20) 

k=l,k£i 

Utilization of the second-order bounds requires evaluation of terms of the form 
P ( Ej ) where E{ denotes the event ’failure in limit state i”. The intersection terms 
refer to domains such as fix shown bounded by the non-linear limit state functions 
g%{X) — 0(i = 1,2,3) in Fig. 3. The individual failure mode probabilities in the 
first-order analysis are determined as 


Pf = <K-0) (21) 

In standardized independent normal X space, the linear limit state function is 
given by 

n 

yi(X) — /?j + E a ij X j (22) 
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where n is the number of random variables. 


The angle between the two limit states provides information about the correlation 
of the two failure modes. The correlation coefficient is obtained as 

n 

Pij = X] OtrOjr = COSt/ij (23) 

r= 1 

Once fii, (3j and are obtained, the computation of joint probability of failure 
can be carried out. Eg. 19 and 20 can be used to compute the second-order bounds 
for system failure probability estimation. 

The above method only provides a tool to approximately estimate the failure 
probability correlation of two different failure modes in a multi-discplinary system. A 
more accurate approach would be the imposition of one failure mode on anther mode 
and reanalysis of the latter. For example, consider two failure modes in a heat transfer 
system: structural failure and heat transfer failure. Structural failure happens when 
the stress, caused by the fluid pressure and temperature difference between outer and 
inner surfaces exceeds limiting value of strength. The heat transfer failure happens 
when the temperature of the contained liquid can not be kept at a certain level. 
When the thermal failure occurs, the increase in the temperature field also causes 
changes in stress field. The structural failure probabilty can be re-estimated under 
this changed stress field and the result can be considered as the interactive failure 
probability under influence of heat transfer failure. A numerical examples will be 
shown for this approach in a later chapter. 
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Symbol 

Meaning of symbol 




Representation of an event 

( 

3 


Representation of an event of a failure 




OR-gate 

1 

i 

n 


AND-gate 


Figure 6: Standard symbols used in fault tree analysis 

Probabilistic fault tree analysis 

NESS US system risk assessment (SRA) uses probabilistic fault tree analysis 
(PFTA). A fault tree is a mathematical construction of assumed component fail- 
ure modes (bottom events) linked in series or parallel leading to a top event, which 
denotes system failure. Standard graphical symbols are used to construct the fault 
tree picture, by describing events and logical connections. These are shown in Fig. 
6, and a simple PFTA is shown in Fig. 7. 

• Fault Trees with a Single AND-GATE 

Consider the fault tree in Fig. 8. Here the top event occurs if and only if all 
the bottom events E\,E 2 , ..., E n occur simultaneously. A system with AND-GATE 
is very similar to a series system structure. 
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Symbols Description 

AND-gate 

_n; — 
Ei 

D 

The AND-gate indicates that the output 
event A occurs only when all the input 
events Ej occur. 

m 


fl 

I 

OR-gate 

i 

The OR-gate indicates that the output 
event A occurs if any of the input events 



II 




[®3 



— 1 h t occur. 

E3 


Figure 8: Fault tree with a single AND-GATE and a single OR-gate 
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Figure 9: Schematic of NESSUS 

• Fault Tree with a Single OR-GATE 

Consider the fault tree in Fig. 8. The top event occurs if at least one of the 
bottom events E\, E ?, ..., E n occurs. The structure of this fault tree is similar to the 
paralell system structure. 

A schematic of Version 6.0 of the NESSUS (Numerical Evaluation of Stochas- 
tic Structures Under Stress) probabilistic structural analysis computer program is 
shown in Fig. 9. As shown in the diagram, the the NESSUS includes other modules, 
namely the System Risk Assessment (SRA) and Simulation Finite Element (SIM- 
FEM) modules. The random field pre-processor (PRE) provides data manipulation 
needed to express the uncertainties in a random field as a set of uncorrelated random 
variables. The user- subroutine which defines the response model (UZFUNC) enables 
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users to define required limit state with the computed response. This study will 
mainly use FEM, PFEM and FPI for reliability analysis. The NESSUS program is 
quite comprehesive with respect to structural reliability estimation. As mentioned in 
Chapter I, the purpose of this study is to develop a technique by which the NESSUS 
program can be used for the system reliability analysis of multi-discplinary systems. 
The following chapters describe this technique in detail. 
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CHAPTER III 


ANALOGY BETWEEN ENGINEERING SYSTEMS 

Introduction 

Since NESSUS/FEM program has been mostly applied only to structural analysis, 
a thermal or a fluid mechanical system needs to be converted through an analogous 
model to a structural system on which the NESSUS program can be applied for 
analysis. Then the probability analysis for a heat transfer system or a fluid mechanics 
system can be carried out by NESSUS. By doing so, a system with heat transfer, 
fluid flow and mechanical stress problems can be analyzed by NESSUS automatically 
with FEM, PFEM, FPI and SRA modules for system reliability analysis. 

In this chapter, a new methodology is presented for one-dimensional steady-state 
heat transfer analysis and one-dimensional steady-state uniform flow problem using 
a structural finite element program. First, the use of the analogous models is intro- 
duced for the analysis of systems involving one-dimensional steady-state heat transfer 
and simple one-dimensional steady-state uniform flow in closed conduit systems. 

Heat transfer analysis through structural analogy 
One-dimensional steady-state heat transfer 

We begin our analysis of one-dimensional, steady-state conduction by discussing 
heat transfer with no internal generation. The objective is to determine the expres- 
sions for temperature distribution and heat transfer rate in common geometries. 
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Figure 10. Differential volume for the derivation of the general equation of heat 
conduction 

The concept of thermal conductivity (analogous to stiffness in stress analysis) is 
introduced as an aid to solving conduction heat transfer problems. Consider a three 
dimensional differential volume shown in Fig. 10, The general heat equation is 




3 t \ . gr 

& > + « = ^ aT 


(24) 


where K is the thermal conductivity of the material. K §J, K §£ are related to 
heat flux in a direction perpendicular to the surface, q is the rate at which energy is 
generated per unit volume of the medium. The density p and specific heat c p are two 
thermodynamic properties. The product pc v is the volumetric heat capacity. pc p ^ 
is the time rate of change of the internal (thermal) energy of the medium per unit 
volume. 
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Figure 11: Heat transfer through a plane wall 

If the heat transfer is one-dimensional and steady state, any differentiation with 
respect to time is equal to zero and there is no internal heat generation, so Eq. 24 
reduces to 


d. T dT . 

li^ KA 3J^ = 0 


(25) 


The heat flux is a constant, independent of x. 

As shown in Fig. 11, a plane wall separates two fluids of different temperatures. 
Heat transfer occurs by convection from the hot fluid at T^i to one surface of the 
wall at T,,i, by conduction through the wall, and by convection from the other surface 
of the wall at T a< 2 to the cold fluid at T ' 00 , 2 . 

Assuming the thermal conductivity of the material to be constant, Eq. 25 may 
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be integrated twice to obtain the general solution 


T(x) — C\X + C 2 


( 26 ) 


To obtain the constants of integration, C\ and C 2 , boundary conditions must be 
introduced. These are: 

T(0) = T til (27 ) 

AL) = T 3i2 (28) 

Applying the condition at x = 0 to the general solution, it follows that 


T,,i — C 2 

Similarly, at x = L 

Tj,2 = C\L + C 2 = C\L + T,'i 

in which case 
A 2 — Ts, 1 „ 

Z 

Substituting into the general solution, the temperature distribution is then 

T «=[". ^ 1 { rl } 

where iVj = 1 - N 2 = £ 

The heat flow can be determined by Fourier’s law, that is 

„ A dT 
q = -Ad— 

ax 


( 29 ) 


( 30 ) 


( 31 ) 


( 32 ) 


or 


q = -KA 


’-(Ti-TtY 


KA 


(Ti - T 2 ) 


( 33 ) 
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Figure 12: Bar subject to tensile force F 
Stress analysis of a bar element 

Now consider a linear-elastic, constant cross-sectional area (prismatic) bar element 
shown in Fig. 12. Using Hooke’s law, the differential equation governing the linear- 
elastic bar behavior is 


sK 


= 0 


(34) 


where U is the axial displacement function in the x direction and S and E are 
cross-sectional area and Young’s modulus of elasticity respectively. 


U = [ Ni N 2 


l h 
U 2 


where Ni = 1 — jj, and N 2 = 

The strain-displacement relationship is 

_ dU_ _ U^-Jh 
£x dx D 


( 35 ) 
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Table 1: Analogous quantities for structural and thermal systems 


Heat transfer 

Structure 

Heat flux q 

Nodal force /i 

Temperature T(x) 

Displacement U(x) 

Inverse of heat transfer resistance 
Conduction: convection: hA 

Structural stiffness ^ 


We obtain 

F = ES ( £l lr L ) <*) 

Also, by the nodal force sign convention of Fig. 12, 

fi = -F (37) 

So Eq. 36 becomes 
E 9 

/, = - U 2 ) (38) 

Analogous modeling between heat transfer and structure 

Comparing Eq. 38 with Eq. 33, the similarities become apparent. These two 

equations indicate a direct analogy between heat transfer and structural analysis. 

The analogous quantities are listed in the Table 1. 

With this analogy, we are able to model a heat transfer problem into a stress 
analysis problem. 
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In the plane wall, we refer heat transfer resistance of conduction R to that 


is 


Rcond. — 


T 3 ,i — T s< 2 


KA 


(39) 


Considering the structural system, Hooke’s law provides stiffness of the form 


, h _ ES 

Ui-U2 D 


(40) 


Comparing Eqns. 39 and 40, and considering ^4 and ^ as analogous qualities, 
fc— d can be traced to be analogous to K. 

A heat transfer factor may also be associated with convection at a surface. From 
Newton’s law of cooling, 


q = hA{T s - T^) 


(41) 


where h is Planck’s constant of convection heat transfer coefficient, T, is the surface 
temperature and Too is the ambient temperature. 

The thermal resistance for convection is then 


R 


conv — 



l 

hA 


(42) 


The equivalent thermal circuit for the plane wall with convection surface condi- 
tions is shown in Fig. 11. The heat transfer rate may be determined from seperate 
consideration of each element in the network, that is, 


q — hiA^Tooj — T Si i ) 


KA 

L 


(T 3A - T,, 2 ) = h 2 A(T aa 


Too, 2) 


(43) 
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In terms of the overall temperature difference, T , - T*,, 2 , and the effective 
thermal resistance R e //, the heat transfer rate may also be expressed 


as 


, _ To C4 - ^00,2 

«.// 


(44) 


Because the conduction and convection resistance are in series and may be summed 
up, it follows that 


R - 1 x 1 ■ 1 

~ I 7 + 77T + 


(45) 


hiA KA h 2 A 

Consider a bar consisting of three different materials which are denoted as ele- 
ments 1,2 and 3. The effective stiffness for this composite bar is 
1 


k e jj — 


JL _l J_ 4 . J_ 
*1 ^2 ^3 


Comparing the above equations Eq. 45 and Eq. 46, the analogy is k eff 
that is 


D 111 
Rejf < > r- + T- + — 

«1 « 2 «3 


Substituting with Eq. 45, we obtain 


1 L 1 
hyA + KA + h 2 A 


1 1_ j_ 

ki k 2 "** k 3 


where 


ki = 


Ex Si 

Dx 


(46) 




(47) 


(48) 


(49) 


fe = 


E2 Si 

d 2 


(50) 


k 3 = 


EzSx 

D* 


( 51 ) 
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Substituting the corresponding terms in Eq.48, we obtain the equivalent quanti- 


ties 


E i 


E 2 



(52) 


(53) 


E 3 <_+ (A a A)(^) (54) 

With these analogous quantities, we use the NESSUS/FEM beam element with 
the E values replaced by the values involving heat transfer problem outlined above. 
The boundary conditions for the bar are the end displacements corresponding to the 
ambient temperature of the wall. After the structural analysis, we get the tempera- 
ture distribution from the corresponding displacement distribution in the output. 

Heat transfer in composite walls 

Equivalence concepts for thermal-structural analysis may also be used for more 
complex systems, such as composite walls and radial heat transfer systems. Fig. 13 
shows a series composite wall. The one-dimensional heat transfer rate for this system 
can be expressed as 


<? = 


Too* ~ 7 * 00,4 

ZR 


(55) 


where To^ — T ^ )4 is the overall temperature difference and the summation includes 
all thermal resistances. Hence 
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Figure 13: Equivalent thermal circuit of a series composite wall 
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Figure 14: Structural analog for the series composite wall heat transfer 


7 * 00,1 ~ Too a 

R\ 4- R'i + Rz + Ri + R$ 


(56) 


7 * 00,1 ~ 7 * 00,4 

(1/M) + ( L a /K a A ) + ( Lb!K b A ) + {Lb/ K B A) + (1 /M) 


(57) 


Alternatively, the heat transfer rate can be related to the temperature difference 
and resistance associated with each element. For example, 


T i rri rp rp rp rjn 

00,1 — * s t l __ ^ s , 1—^2 __ ^2 — ^3 

q ~ (1 /h,A) “ (La/KaA) ~ (LbIKbA) 


(58) 


The analogous structural model for this series composite wall heat transfer prob- 
lem is shown in Fig. 14. The bar consists of five elements with stiffnesses of 
ki,fc 2 ,k 3 ,k 4 ,k 5 . Using the mechanical structure equivalence for convection and con- 
duction, we obtain 


E={ 


(M) (?) 

(¥) (?) 


for convection 
for conduction 
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Engine oil 



/>l Q> 2<> 3 Q) 405 Q i 6 


(b) Structural model 

Figure 15: A heat exchanger for engine oil and refrigerant fluid 
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Heat transfer in radial systems 


Cylindrical and spherical systems often experience temperature gradients in the 
radial direction only and may therefore be treated as one dimensional. 

Fig. 15 shows an example of a heat exchanger, whose inner cylinder is used to 
store engine oil and the outer cylinder is used to transfer the refrigerant fluid to 
cool down the oil temperature. The outer insulated covering is assumed to isolate 
the system from the ambient environment. For steady state conditions with no heat 
generation, the appropriate form of the heat equation is 


TTr (*'£) = 


0 


(59) 


The rate at which energy is conducted across any cylindrical surface in the solid 
may be expressed as 


„ dT r ,dT 

" = ~ KA i; = -K^rL)- 


(60) 


where A — 2nrL is the area normal to the direction of heat transfer. 


The thermal resistance is 


R'fJ = 


1 


ln(r 2 /r x ) 1 

l r* TS T ' 


hi2irriL 2nKL 


(61) 


which includes both conduction and convection. 


The heat transfer rate for a unit length of the cylinder therefore is 


Q = 


Too,! — Too ( 2 


i 


+ 


ln(ra/n) 


h\2rcr\L ' 2nKL 


+ 


(62) 


/l2 27Tf2 X/ 
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Table 2: Analogous quantities for heat transfer in a radial system. 


Structure 

Heat transfer 

Ei 

( 2?rr iM (§■) 


e 2 



e 3 

(271-7-2*2) (ff) 



The structural analog for the cylinder is shown in Fig. 15(b). In this case, the 
mechanical structure equivalence is k e fj 4 — >• jQ 1 • The corresponding equivalent 
quantities are listed in Table 2. The E values are input to NESSUS structural anal- 
ysis, and the output displacements from NESSUS give the temperature distribution. 

Numerical example for heat transfer solved with NESSUS/FEM 

Fig. 17 shows the sectional view of the cylindrical copper heat exchanger which 
the engine oil flows through. The copper wall thickness is 0.281 in. The radius 
to the surface of the insulation pipe covering (fc; = 0A28Btu/(h - ft 2 F)) is 
1.33 in. The fluid in the outer container is controlled at a constant temperature 
of 70° F. The forced convection heat transfer occurs between the outer surface of 
the insulation covering and the flowing fluid with h = 10 Btu/(h - ft 2 F). The 

surface temperature at the insulation covering is 35.298° F. The structural analogy 
model is used to determine the inside temperature of the tube, assuming steady 
state, one dimensional, uniform properties in each material, forced convection cooling 
and negligible thermal radiation. The conductivity coefficient of copper at room 
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' ? 



T F - 70°F R kI 7i /?* 2 r 2 /? c 73 = 57 

(a) Thermal circuit 


k3=E 3 S 3 /D 3 ki = E. S j/Dj 



(b) Structural analogous model 


Figure 16: Analogous model for the heat exchange 

temperature is K c = 223 Btu/(h — ft 2 —° F). 

The thermal circuit is shown in Fig. 16(a). Fig. 16(b) shows the structural 
analog model for this heat transfer problem. Beam element type 98 in NESSUS/FEM 
element type library is adopted. Three elements represent three heat transfer forms 
involved in this problem, which are the forced convection between the surface of 
the insulation covering and the ambient air, conduction through the copper layer, 
and conduction through insulation covering, respectively. Therefore, in terms of the 
structural model, we must assign three different material elastic constants for this 
beam structure. Since the NESSUS/FEM utilizes the Nodal-based data input , two 
duplicate nodes are used at each boundary between elements 1 and 2, and between 
elements 2 and 3. The room temperature 70 °F becomes the boundary displacement 
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r 0 = 1.330 in 
T) = 0.669 in 
r 2 = 0.950 in 


7» = 70° F 

Ji = 10 Btu/(h - ft -°F) 
k = 0.428 Btu/(h- ft 2 -°F) 


Figure 17: Section of a radial heat exchanger 

70.0 at node 1. A concentrated load F at point 3 is -241.66 lb. The length of the 
structural element is 10.0 units, the sectional area is 1.0. It should be noted that the 
units used in the structural model here do not have any real meaning in terms of a 
real structure. They are simply used to facilitate the structural analysis. 

The equivalent values are calculated as follows. 


. ( 2t tK c \ 


2 7r x 223 

w in 

: 39955.476 

(63) 

U n ( r 2/ r l)J 

‘ \Si) 

ln(0.95/0.669) 

. ( 2rcK i \ 


2tt x 0.428 

79.923 

(64) 

Vln(r 0 /r 2 )J 

1 \sj 

ln(1.33/0.95) x 10 ~ 

: (2irr 0 h) 


< 11.33/12 x 10.0 x 10 

= 69.63867 

(65) 

2tt roA(Too - 

T 3 ) = 27T 

x 1.33/12 x 10(35.298 

- 70) = -241.66 

(66) 
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NESSUS/FEM uses this data, and gives the output of the displacement distribu- 
tion in the structure as follows: 

U x = 70.000 
U 2 = 35.298 
U 3 = 35.298 
U 4 = 5.0614 
U 5 = 5.0614 
U 6 --- 5.0009 

Converting the above displacement information to the equivalent temperature 
distribution, we obtain: 

^ = 70.000° F 
T x = 35.298° F 
T 2 = 5.0614° F 
T 3 = 5.0009° F 

The data and the output files are shown in Appendix A. 

Fluid flow analysis through structural analogy 

Equation of motion for fluid flow 

The Bernoulli equation gives a relationship between pressure, velocity, and posi- 
tion or elevation in a flow field. Normally, these properties vary considerably in the 
flow, and the relationship between them if written in differential form is quite com- 
plex. The equation can be solved exactly only under very special conditions. There- 
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L 


Figure 18: Fluid in a constant diameter duct 
fore, in most practical problems, it is often more convenient to make assumptions to 
simplify the descriptive equations. The Bernoulli equation for steady, incompressible 
flow along a streamline with no friction (no viscous effects) is written as [10] 
p V 2 

- + T + gz = C (67) 

where 

p is fluid pressure 
p is the desity of the fluid 
V is the flow velocity 
g = 32.174 ft/s 2 , and z is height. 

For a horizontal pipe shown in Fig. 18, Zt = z 2 . From continuity, A t Vi = A 2 V 2 . 

Because D x = D 2 , then Ai = A 2 , and therefore, V\ = V 2 . The Bernoulli’s equation 
reduces to 

Pi=P2 (68) 
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1 w Pdz 

■<= — 


pA 


(p+dp)A 




Figure 19: Control volume of a system: flow in a duct. 

This result is not a proper description of the situation, however. For flow to be 
maintained in the direction indicated in Fig. 18, pi must be greater than p-i in an 
amount sufficient to overcome friction between the fluid and the pipe wall. In order 
to apply Bernoulli’s equation and obtain an accurate description, we must modify 
the equation with a friction term. 

Consider flow in a pipe as shown in Fig. 19. A control volume that extends to 
the wall (where the friction force acts) is selected for analysis. 

Note that a circular cross section is illustrated, but the results are general until 
we substitute specific equations for the geometry of the cross section. The forces 
acting on the control volume are pressure normal to the surface and shear stress 
acting at the wall. The momentum equation is [10] 

= f Jv zP V n dA (69) 
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where 

K is the fluid velocity along the longitudinl direction 
V n is the norminal fluid velocity 

Since the flow out of the control volume equals the flow in, the right-hand side 
of this equation is zero. The sum of the forces is 


pA - r w Pdz - (p + dp)A = 0 


(70) 


where 


A 

Pdz = 


= cross-sectional area 

the surface area (perimeter times length) over which the 
wall shear t w acts 


The equation reduces to 


r w Pdz + Adp = 0 

Rearranging and solving for pressure drop, we get 
dp _ At^ 

dz D h (72) 

We have thus expressed the pressure drop per unit length of the conduit in terms 
of the wall shear and the hydraulic diameter. Eq. 72 is a general expression for any 
cross section. It is convenient to introduce a friction factor /, which is customarily 
defined as the ratio of friction forces to inertia forces: 

, _ 4tu, 

1 \pV a (73) 
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Figure 20: Laminar flow in an annulus 


where V is the average flow velocity. 

By subsitution into Eq. 72, we obtain 


dp = 


pV ' 2 fdz 

2 D h 


(74) 


Integrating this expression from point 1 to point 2 a distance L apart in the 
conduit yields 


V 2 = --jyAp (75) 

pfL 

Eq. 75 gives the relationship between the velocity and the pressure drop in the 
duct due to friction. This equation can be applied to two flow regimes - laminar and 
turbulent flow. However, caution must be excercised when determining the friction 
factor /. 

This equation can also be applied to flow through noncircular cross section such 
as rectangular duct and annulut Fig. 20 shows the laminar flow in an annulus. 
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Table 3: Analogous quantities between structural and flow systems 



The annulus flow area is bounded by the inside surface of the outer duct (radius 
Ri) and the outside surface of the inner duct (R 2 ). We define the ratio of these 
diameters as 



(76) 


in which 0 < k < 1. 

The friction factor used in Eq. 75 is defined as [10] 

1 _ R e [1 +k 2 l + k' 

f 64 [ 1 k + ln(fc)J ( 77 ) 

where f? e = ^1(1 - fc). 

Compare Eq. 75 with Eq. 38 concerning the beam structure subjected to the 
end nodal force, as discussed in previous section: 

E ^ 

h = 2 ) (78) 


We are now able to set up the analogous qualities listed in Table 3. 
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Numerical example of flow in a tube solved with NESSUS/FEM 


Consider the refrigerant flow in a copper tube as an example to demonstrate how 
NESSUS/FEM can be applied to problems in fluid mechanics. 

A horizontal copper duct as shown in Fig. 18 with inside radius of 0.669 in, and 
1,200 in in length. If the inflow pressure pi is 1838.7 psi, assuming the refrigerant 
is Freon F-12 under a temperature of 5 °F, p is 0.0499 lb/ in 3 . The friction factor / 
is assumed to be 0.03, and V is 15.5 in /sec. The objective is to obtain the outflow 
pressure p 2 using NESSUS/FEM. 

First of all, we need to identify all the equivalent quantities for structural analysis. 
We assume a single element beam structure subjected to a concentrated force equal 
to -240.25 units. The beam element has a section of 0.1 mx 0.1, and a length of 1.0. 
The boundary condition is an initial displacement of 1838.3 units at node 1. Again, 
it should be mentioned that the units used here do not have real meaning in terms 
of a real structure. According to Table 3, the analogous quantities can be obtained 
as follows 



2 x 2 x 0.669 
0.0499 x 0.03 x 1200 


x 10 = 14.896 


/i = V 2 = 15.5 2 = 240.25 


F = -240.25 


This data is input to NESSUS/FEM, and the displacement at point 2 is ob- 
tained as 1677.4. Converting this displacement to the fluid model, we get the output 
pressure p 2 = 1677.4 psi. 

The data and the output files are shown in Appendix B. 
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CHAPTER IV 


MULTI-DISCIPLINARY SYSTEM RELIABILITY ANALYSIS 

Introduction 

After transforming the heat transfer and fluid mechanics problems into corre- 
sponding structural analog models and using NESSUS/FEM to perform the finite 
element analysis, we can define the individual failure modes in NESSUS/FPI. Then 
NESSUS/PFEM can be employed to integrate FEM and FPI programs to obtain 
the failure probability and CDF for each failure mode. The failure mode for heat 
transfer problem would be defined as, for example, the event that the temperature 
at a certain location is lower or higher than the required temperature. The failure 
mode for fluid flow would be defined as the flow pressure exceeding a certain pressure 
level, and the structural failure is defined as the stress exceeding either the ultimate 
strength or the yield strength of the material. 

Upon the completion of failure probability analyses of individual failure modes, 
the system failure analysis can be pursued. The different failure modes involved in 
a system have different impacts on the overall performance of a system. Some types 
of failure such as structural failure are critical to the system. If the material used to 
construct the main parts of the system fails, the whole system can no longer function. 
Such failure is called critical failure. Other failures modes such as thermal failure of a 
heat exchager do not destroy the system bu’ degrade the performan x of the system. 
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Such failure is referred as to functional failure. The function of fluid flow will fail 
when the outflow pressure rises higher than the designed value, but the system can 
still be working until the pressure increases to the level which will cause the system 
to shut down. The individual failure modes can also be correlated to each other. For 
example, the temperature field in the thermal failure mode affects the stress field in 
the structural mode. The flow pressure definitely has impact on the stress. However, 
in some cases, the component-level events in a system is considered as independent 
events. In the example which will be discussed later in this chapter, the thermal 
failure and fluid flow failure modes do not share correlated input parameters, so they 
are considered as independent of each other. 

Using the analogy method, the thermal and fluid flow problems are analyzed sim- 
ilar to the structural model by means of NESSUS. For physically correlated events, 
the failure mode of one event is imposed on the other. In a system consisting of 
structural, thermal and fluid flow modes, the thermal and flow failures are imposed 
into the structural failure analysis to study the impact of correlated events. The 
failure probability of the whole system is then estimated based on the output from 
the above analyses. 

Individual failure analysis 

Structural failure mode 

For a copper duct of a heat exchanger shown in Fig. 17 in Chapter III, the 
structural failure mode is defined as that when the tensile stress exceeds the yield 
strength /j,, e /d = 8.0 ksi. In finite element modeling, we use Element type 153 in the 
NESSUS/FEM file. This element is a four-noded quadrilateral lying in the global 
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zr — plane which is defined by cylindrical coordinates. 

This structure is subjected to two types of load - fluid pressure from the inside 
flow and the stress caused by temperature difference between the outer surface and 
the inner surface. 

One convenient feature of the NESSUS/PFEM is that we can impose different 
temperatures at the inner surface of the pipe and obtain the different probability 
results under different temperature conditions. This enables us to investigate the 
effect of different temperature levels on the structural failure probability. By doing 
so, the relationship between the failure modes in two dicsiplines - heat transfer and 
structural mechanics - is established. This is a significant step toward the system re- 
liability analysis with physically correlated failure modes. This will be demonstrated 
in a later section. 

For this structural failure model, we first suppose that the temperature failure 
(which will be descibed in the following section) did not occur, that is, the tempera- 
ture at the inner surface of the duct is below 5.0009° F. Given that the inner surface 
temperature is 4.0000°F, using the FEM we obtain the outer surface temperature as 
4.0605°F. 

Also, we assume that the outlet flow pressure is under 1677.4 psi which enables 
the system to work properly. We assume that the outlet flow pressure is 1577.4 psi. 

We input this temperature and flow pressure profile in the structure FEM data 
file, with the random variables defined in Table 4. 
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Table 4: Random variables for structural model 


Variable 

Mean value 

Distribution 

C.O.V. 

Pressure 

P 



0.2 

Modulus of elasticity 
E 

1.7xl0 7 psi 

Normal 

0.1 

Coefficient of thermal 
expansion a 

9.5 xKT 6 

Normal 

0.1 


Table 5: CDF corresponding to different tensile strength levels 


Z-level 

(strength) 

CDF 


mmm 

-206.89941 psi 

0.00000017 

4385.4964 psi 

0.81593991 

558.49989 psi 

0.00002067 

5150.8957 psi 

0.97128351 

1323.8992 psi 

0.00096767 

5916.2950 psi 

0.99813412 

2089.2985 psi 

0.01786435 

6681.6943 psi 

0.99995188 

2854.6978 psi 
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The result attached in Appendix C indicates that the structural reliability when 
the heat exchanger is working properly in thermal and fluid aspects is 0.99999999. 
The failure probability is expressed as 1 — Preliability Therefore, the structural failure 
probability is 1.0 x 10 -9 . 

The key word response type in the NESSUS input data file, FPI section, is set 
equal to 3 which means that the response quantity used in limit state function is 
stress. The corresponding keyword analysis type in FPI section is first set equal to 
1 which means that the probability analysis is for a single Z-level. The Z-level in this 
case is 8,000 psi. The probability result will be under the condition of a < 8,000 
psi, i.e., the structural reliability of the system under certain thermal and fluid flow 
working conditions. 

The CDF is obtained by using PFEM by setting analysis type in FPI section 
equal to 0 which automatically generates a set of different values of Zq (i.e., Z levels 
for a series of stress valus) for probability analysis. The CDF values corresponding 
to different strength Z-levels are shown in Table 5. 

The CDF chart is shown in Fig. 21. It should be noted that the first line of 
the data which contains negative Z-level is eliminated because negative stress is 
considered impractical in this model. The input and output files are attached in 
Appendix C as well. 
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8000.0 


Figure 21: CDF of structural reliability of refrigerant duct 


Table 6: Random variables for thermal model 



Mean value 

Distribution 

Coefficient of variation 


223.0 

Normal 

0.1 

K, 

0.428 

Normal 

0.1 


11.3 

Normal 

0.1 
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CDF 
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Failure mode in heat transfer 


Using the same example of a heat exchanger as in Fig. 17, we define a failure 
event when the inside temperature is higher than 5°F, because the refrigerant will 
not function properly beyond 5 °F which is considered as a failure in the device we 
studied. First the input data for NESSUS/PFEM is set up to obtain the reliability 
under this failure mode, then the data file is set up with different Z - levels to obtain 
the CDF, which provides reliability estimatecorresponding to different temperature 

levels. The random variables K c and K, and h for heat transfer are defined in Table 

6 . 

In order to use NESSUS/PFEM, the analogous quantities E\, E 2 , E 3 and F 
are calculated from Eqs. 63, 64, 65, 66. Because the distribution of the random 
variables K c and A, and h is normal and Ej, E 2 , E 3 are linear to K c and A", and 
h, the distribution of random variables E u E 2 , E 3 is also normal. The mean values 
and standard deviations of of E u E 2 and E 3 are input to NESSUS/PFEM. 

The Z-level is 5.0°F, so P(Z < Z 0 ) is the probability the device can keep the 
inside fluid temperature under 5.0°F, which is the thermal reliability of the system. 
We set up the keyword in FPI section analysis type equal to 1 which means the 
probability analysis is performed for a single Z-level. The result is attached in Ap- 
pendix D. The thermal reliability of this device is 0.9099214. Therefore, the thermal 
failure probability is 9.00786 x 10~ 2 . 

The CDF is obtained by setting up the FPI keyword analysis type equal to 0 
which automatically generates a set of different values for Z 0 (i.e., Z-levels). The 
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Table 7: Random variables for flow model 



Mean value 

Distribution 

Coefficient of variation 

D h 

0.669 

Normal 

0.1 

V 

18.71 

Normal 

0.1 


CDF output is shown in Appendix D and the CDF curve is shown in Fig. 22. 

Failure mode in fluid flow 


Next we consider the one-dimensional fluid flow in a duct of a heat exchanger. 
The failure mode is defined as the pressure at a certain point along the duct rising 
above the value at which the system cannot function properly. 

The example of a duct in a heat exchanger shown in Fig. 16 is used. The only 
difference is that V is assumed to be 18.7 in j sec. We define the failure mode when 
pressure rises above 1677.4 psi. The Z-level is therefore 1677.4 psi. The keyword 
response type is set as 1 for the displacement output which is the analogy of the 
pressure. The random variables related to fluid flow are defined in Table 7. The 
analogous quantities for use in NESSUS/PFEM are calculated according to Table 3 
as follows: 



2 x 2 x 0.669 
0.0499 x 0.03 x 1200 


x 10 = 14.896 


/i = V 2 = 18.7 2 = 350.0 


F = ~f x - -350.0 
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CDF 




Since £ is a linear fuction of Dh, its distribution is also normal. But since /i = V 2 , 
the distribution of f\ is actually chi-square (x 2 )- However, we have used the normal 
distribution for f\ in this study as an approximation. The friction fator / is assume 
to be a constant. 

A NESSUS/PFEM input data file is compiled. The reliability is obtained by 
setting the keyword analysis type equal to 1, and the CDF is obtained by setting it 
to a value of 0 which automatically generates a set of different values for Zq. The 
PfEM input and output files are shown in Appendix E and the CDF curve is shown 
in Fig. 23. 

The reliability is 0.98670241 and therefore the failure probability of the output 
flow pressure being higher than 1677.4 psi is 1.329759 x 10 -2 . 

Multi-disciplinary system reliability 

After the individual failure modes are identified and analyzed, the system reli- 
ability analysis can be pursued. Fig. 24 shows a device which is used to transfer 
refrigerant fluid through a copper duct. The duct is installed in an enclosed cham- 
ber which is maintained at a constant temperature of 70 ° F. The thickness of the 
copper wall is 0.281 in. The radius to the surface of the insulation pipe covering ( 
mean value of K{ equals to 0.428 Btu/(h — ft 2 — 0 F), c.o.v. equals to 0.1 ) is 1.33 
in. Forced convection heat transfer occurs with h — 10 Btu/(h — ft 2 — 0 F) ( mean 
value with c.o.v. equals to 0.1). The thermal conductivity of copper K c has a mean 
value of 233.0 Btu/(h — ft 2 — ° F) with a c.o.v. 0.1. The surface temperature at 
the insulation covering is 35.298° F. The inflow pressure p\ is designed to be 1838.7 
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psi. The refrigerant is Dichlorodifluoromethane ( Freon F-12 ) under a temperature 
of 4°F. p is 0.0499 lb/m 3 . The friction factor / is assumed to be 0.03, and V is 
18.71 in j sec ( assuming V 2 is normal distribution and has a c.o.v. of 0.1). The inner 
radius of the duct has a mean value of 0.669 in and a c.o.v. of 0.1. 

The above data are the same as the model shown in Fig. 17 of Chapter III, which 
are used in FEM analysis and reliability estimation of individual failure modes. The 
system shown in Fig. 24 is simply the combination of the previous individual models 
which have been analyzed in different disciplines. The following is to demonstrate 

how the analysis results of the individual failure modes can be integrated into the 
analysis of a whole system. 

The failure of the system consists of the individual failure modes in three disci- 
plines: structural failure, thermal failure and fluid mechanics failure. 

• First of all, the duct should work without any damage to the structure, i.e. the 
duct should be structurally sound without yied or crack. If yield occurs, then 
structural failure is assumed to occur. We denote the structural failure as E\. 
The structural failure is a critical failure in this system. 

• The refrigerant liquid this device transports is sensitive to temperature changes. 
The requirement is that the temperature cannot be higher than 5°F for the 
next process to proceed. If the temperature of the liquid rises higher than 5° F, 
then thermal failure occurs, which we refer to as E 2 . The thermal failure is a 
non-critical functional failure in this system. 

• It is required that the fknd flow be maintained a certain presure at the 
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Input valve 

Outlet to the next process 




Figure 24: Refrigerant duct through a chamber 
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ends which enables the refrigerant fluid to maintain a steady speed to provide 
constan volume in the device. If the flow pressure rises higher than 1677.4 
psi, the flow failure occurs, which is referred to as E 3 . The flow failure is a 
non-critical functional failure in this system. 

From the previous section, the probability of the individual failure modes of the 
system shown in Fig. 24 have been obtained, wich are: 

P(Ei) = 1.0000000 x 10 -9 
P(E 2 ) = 9.0078600 x 10 -2 
P(E 3 ) = 1.3297590 x 10" 2 

As was indicated in the previous section, we can also impose one failure mode 
upon the other. In this case we can impose the thermal failure ( which happens 
when the fluid temperature rises above 5 ° F ) and the fluid pressure failure ( which 
happens when the outlet fluid pressure rises above 1677.4 psi) upon the structure 
respectively. In the FEM file for the structure model, the corresponding data are 
modified to impose those failures. 

First, we assume that the temperature failure occurs while the fluid pressure is 
still lower than 1677.4 psi , say 1577.4 psi, i.e., the fluid flow is operation in safe mode. 
By redefining the temperature profile in the FEM data deck as 6.06°F (for example) 
in the inner layer of the wall and 6.12°F in the outer layer of the wall, which mean the 
thermal failure occurs, we inpose the thermal failure to the structural model. The 
PFEM result gives us the structural failure probability under the condition that the 
thermal failure occurs. In this case, the structural reliability is 0.99999996, therefore. 
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the failure probability P{E\I E 2 ) = 4.0 x 10~ 8 . 

Now we impose the fluid pressure failure upon the structural model. The fluid 
mechanical failure occurs when the pressure at the outlet rises above 1677.4 psi, say 
1777.4 psi. The structural FEM file is modified by redefining the pressure profile 
according to this failure pressure. It should be noted that the temperature profile 
should remain under the normal working condition, which is that the temperature 
in the inner layer of the wall of duct is under 5°F, say 4°F, i.e., the thermal as- 
pect of the system is operating in the safe zone. The result indicates that in this 
case, the structural reliability is 0.99999633, therefore, the conditional probability, 
P(Ex/Es ) = 3.67 x 10- 6 . 

Next, both thermal and fluid mechanical failures are imposed that is, the fluid 
temperature rises above 5 °F, and the outlet flow pressure rises above 1677.4 psi. 
Modifying the input FEM data deck in structural PFEM file with inner surface 
temperature of 6.06 °F, and the fluid pressure of 1777.4 psi, we can get the result of 
the structural reliability of 0.99999227, which means, P(Ei/ E 2 E 3 ) — 7.73 x 10 -6 . 

The conditional probabilities of structural failure have been obtained as 

P(Ei/Et) = 4.00 x lO" 8 

P(Ei/E 3 ) = 3.67 x 10“ 6 
P(E x IE 2 E 3 ) = 7.73 x 10~ 6 

System reliability computation 

System reliability analysis can be performed in two different ways, depending on 
the definition of the system failure. In the first (traditional) method, we define that 
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System Failure 



Structural Thermal Fluid flow 
Failure Failure Failure 


Figure 25: Fault tree for the system with three critical failures 
the system failure occurs when any component-level failure occurs. In the system 
mvoving structural, thermal, fluid mechanical failure modes, i.e., E u E 2 and E 3 , the 
system failure can be illistrated in a fault tree shown in Fig. 25. 

As discussed in Chapter II, the probability of system failure P(E) can be obtained 
using the following equations: 

P(E) = P(E\ U£ 2 U E$) 

The above expression can be expaned as: 

P(E) = P(E r ) + P{E 2 ) - P(E\ n E 2 ) + P(E 3 ) - P{Ex D E 3 ) 

-p{E 2 n e 3 ) + P(Ei r\E 2 n e 3 ) (80) 

Since the joint probability is not always available, an approximate method is to 
consider the individual failure modes as independent and ignore the correlations. In 
our case however, the conditional probabilities have been calculated. Therefore, the 
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system failure probability can be computed as 


P(E) = P(E l ) + P(E 2 ) + P(E 3 ) - P(E l /E 2 )P(E 2 ) - P(E 1 /E 3 )P(E 3 ) 

-P(E 2 /E 3 )P(E 3 ) + P{E x IE 2 E 3 )P{E 2 E 3 ) (81) 

For systems involving many failure modes, approximation methods are used to 
predict the system failure probability (or reliability), such as first-order bounds or 
second-order bounds [1]. 

Since no correlation is assumed between the failures of flow pressure and fluid 
temperature, we assume P(E 2 fl E 3 ) is equal to P(E 2 )P(E 3 ) in our analysis. 

Substituting the numerical results from the previous discussion into Eq. 81, 
we obtain the probability of the system failure of the heat exchanger, P(E) = 
0.10217832. 

As mentioned before, the above failure probability is an estimation of system 
failure in case any failure occures which includes both critical and non-critical func- 
tional failures. Now we will pursue the probability estimation for the system critical 
failure which, in our case, is structural failure. During the service cycles, the thermal 
and fluid mechanical failures may be non-critical, i.e., their occurence does not cause 
total system failure. They will cause the system to fail in some functions as designed, 
such as keeping the fluid under certain temperature or keeping the outlet fluid pres- 
sure under certain value. However, if the the system keeps operating, the changes in 
temperature and fluid pressure will cause progressive damage to the structure due 
to load redistribution. The estimation of critical structural failure of the system has 
to consider the progressive damage caused by all components. Structural reanalysis 
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Table 8: Structural failure probability under various temperatures 


T(°) 

0.0 

5.0 

10.0 

15.0 

P(Ei) 

0.0 


m 

0.0 


is used to account for the effect of non-critical damage on critical failure mode. In 
the refrigerant model, a reanalysis procedure is performed to accurately estimate the 
failure region segments for structural failure mode affected by progressive damage 
caused by thermal and fluid pressure changes within the system. The overall struc- 
tural failure probability is obtained through the union of the failure region segments 
defined by each limit-state function. 

We can also impose a series of temperatures under which the system may be 
operating upon the structural model to examine the temperature impact on the 
structural failure probability. Just as we did before, the failure probability is obtained 
as (1 - Reliability). In this case, we still assume inner fluid pressure is 1577.4 psi, 
which means that the fluid flow mode of the system is operating in the safe zone. 
The results are shown in Table 8. 

Table 8 shows that when the fluid pressure is not considered as a random variable 
in perturbation for probability analysis, the temperature changes do not have a 
significant impact on structural reliability of the system. 

We can also get the structural failure probability under different pressure condi- 
tions by defining a series of the pressure profiles in the FEM data deck for structural 
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Table 9: Structural failure probability under various pressures 


Pressure 

(psi) 

2700.0 

2900.0 



3300.0 

3400.0 

3420.0 

P(E,) 

0.000 

1.000 

xlO -9 

1.920 

xlO -6 

7.687 

xlO -3 

9.195684 

xl0~ 2 

4.06468 

xlO" 1 

4.92622 

xlO" 1 


model. The results are listed on Table 9. In this case, we assume that the inner tem- 
perature is 4°F. 

It should be noted that the NESSUS/PFEM input file gfun.dat is different from 
the previous structural PFEM file in which the pressure is defined as a random 
variable. In the program [Ang and Tang, 1984] to calculate the union of the re- 
gion segments, the random variables once defined can not be changed for different 
limit states. Since in the pressure profile the different pressure levels are presented, 
pressure should not be defined as a random variable. Therefore only two random 
variables are involved in gfun.dat - modulus of elastisity E and coefficient of thermal 
expansion a. The gfun.dat and various gfun.mov files are shown in Appendix F. 

There are two ways of quantifying the effect of progressive damage on critical 
failure. The first is simply to compute the variation of critical failure probability 
with respect to progressive damage. This is shown in Fig. 26 for various pressure 
levels. 

An alternate way is to compute the progressive damage on overall critical failure 
probability. If each critical failure limit sate segment for each progressive damage 
gives the event E,-, then the overall critical failure probability is P(U£_ x Ej), where n 
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Failure probability 
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is the number of limit states ( n = 7 in this case ). 

The output files gfun.mov provide g-functions defining different limit-states for 
each level of damage as 

<K = # + X>,< (82) 

t=i 

From Table 8, it is clear that temperature variations do not have any significient 
effect on the structural failure probability. Therefore, only pressure variations are 
considered as follows: 

The parameters provided by gfun.mov for the structural limit states correspond- 
ing to fluid pressure profile are as follows: 


Pressure 

0 

<*i 

C*2 

2700.0 psi 

6.825696 

0.999887 

-0.015057 

2900.0 psi 

5.721588 

0.999886 

-0.015088 

3000.0 psi 

4.619647 

0.999886 

-0.015071 

3200.0 psi 

2.423429 

0.999887 

-0.015038 

3300.0 psi 

1.328801 

0.999887 

-0.015021 

3400.0 psi 

0.236641 

0.999888 

-0.014959 

3420.0 psi 

0.018494 

0.999887 

-0.015003 


The above data provides parameters for 7 g-functions. Using the above data 
to calculate the union of the region defined by a group of g-functions [Ang and 
Tang, 1984], we obtain the probability of structural failure involving the progressive 
damage caused by fluid pressure. The probability is defined by a lower and upper 
bounds, which in this case are both 0.49261 1 . The data and output files are shown in 
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Appendix G. Those bounds provide the overall structural failure estimation when the 
system experiences various levels of progressive damage. The fluid pressure changes 
in a range from 2700.0 psi to 3420.0 psi. 

The above two methods provide practical tools for multi-disciplinary system re- 
liability estimation using NESSUS. With multiple impositions of one mode on the 
other mode, a close approximation to the failure domain can be constructed, and the 
critical failure probability can be obtained through the union of the failure region 
defined by the various limit-states. 
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CHAPTER V 


CONCLUSIONS AND RECOMMENDATIONS 
Conclusion 

This report has demonstrated the application of equivalence concepts to the re- 
liability analysis of multi-disciplinary systems using NESSUS. A thermal-structural- 
fluid system is used to illustrate the proposed methodology. The analogous model 
is a very powerful tool to analyze the one-dimensional steady state problem in heat 
transfer and fluid mechanics by converting those models into a structural model. 
Then the NESSUS probability analysis program can be implemented and the precise 
system reliability can be evaluated. Both traditional and progressive system fail- 
ure probability methods using NESSUS provide practical tools for multi-disciplinary 
system reliability analysis. 

Recommendations for future research 

This research project demonstrated how the NESSUS program could be applied 
for reliability analysis of engineering systems involving different disciplines, such as 
structure, heat transfer and fluid mechanics. The current models are based on the 
condition of one-dimensional, steady-state for both heat transfer and fluid mechanics. 
More complex systems could be treated in the similar way. However, the scope of 
application of this methodology is largely dependent on the ability of NESSUS/FEM 
to deal with problems in different disciplines under more complicated situations, for 
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instance, a thermal or a flow model in two or three dimensions and non-steady state 
conditions. For more sophisticated systems, the need for separate FEM program 
may be inevitable. Either a new source code should be developed or the existing 
commercial softwares could be integrated into the program. Several users of NESSUS 
have already integrated its FPI module to other FEM analysis program such as 
ANSYS and NASTRAN. Nevertheless, the use of equivalent concepts helps to obtain 
a quick estimate of multi-disciplinary system reliability through the use of NESSUS. 

The idea of progressive damage imposition to quantify nonlinear system reliability 
effects has previously been applied to structural mechanics problems [1]. This study 
extends this concept to multi-disciplinary systems. This appears to be a practical 
methodology for system reliability analysis when failure modes (even in different 
disciplines) have physical relationship with each other. The methodology should be 
pursued further for application to other, more complicated engineering systems. 
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Appendix A 


FEM FILES FOR HEAT TRANSFER OF A HEAT EXCHANGER 


(Refer to Chapter III, Fig. 16) 
NESSUS/FEM input file 


•xp2fem.dat 


Wed May 14 20*37*04 1997 


, HEAT TRANSFER TOR A REGRIOERANT PIPE 


PARAMETER DATA 


20 


•BOITN 
•CONS 0 

•DISP 

•NODS 6 

•DUPL 4 

• BLEU 3 

98 

•forc l 

•PERT 4 4 

•PAIN 

•BEAM 

•END 

c .... model DATA 
•COOA 


3 

4 
3 
6 

•ELEM 

1 

2 

3 

•DOPL 

3 

5 

•BO UN 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0000 

10.0000 

10.0000 

20.0000 

20.0000 

30.0000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


1.0000 

1.0000 

1,0000 

1.0000 

1.0000 

1.0000 


6 

1 

0.0000 





« 

3 

0.0000 





1 

1 

0.0000 





1 

2 

70.0000 





1 

3 

0.0000 





1 

4 

0.0000 





1 

5 

0.0000 





1 

6 

0.0000 





•BEAM 

1 






1 

€ 

1,0000 

1.0000 




•ITERA 

0 

4 





40 

0.0500 





•PROP 

1 

98 

2 

1.0000 

69.63867 

0.3000 

0.0000 

2 . 5E-4 

3 

4 

1.0000 

79.9230 

0.3000 

0.0000 

2.5E-4 

5 

6 

1.0000 

39955.5476 

0.3000 

0.0000 

2.SE-4 

•FORC 







6 

2 

-241 . 64 





•PR1N 








TOTA 

NODE 
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NESSUS/FEM output file 


exp2fe».out Wed May 14 20:50:14 1997 


NN NN 

NNN NN 

NN N NN 

NN N NN 
NN N NN 
NN NNN 

NN NN 


EEEBBEEB SSSSSSS 
*E 55 

« SS 

eSEEBE SSSSSS 
SS 

EE SS 

BCEEEEEE SSSSSSS 


SSSSSSS uu uu 

SS uu uu 

SS uu uu 

ssssss uu uu 

SS uu uu 

SS uu uu 

sssssss uuwuu 


SSSSSSS 

SS 


ssssss 

SS 


NNNNN NNHN 

N N 

NN NNN NN 
NNNNNNN 
NNN NN NN 
N N N 

N N N 

N N N 


FFFFFPFF 

FF 

FF 

FFFFFF 

FF 

FF 

FF 


EEEEEEEE MM MM 
EE MKM MMM 
EE MM H M MM 
EEEEEE MM M MM 
EE MM KM 
BE MM KM 
EEEEEEEE MM MM 


DATE: 14- 5-1997 20:50 - LEVEL 6.0 (109) 


DATED June, 30 


C . . . MEAT TRANSFER FOR A REORIOERANT PIPE 


1992 


•* PERTURBATION NO. 4 ** 


SUMMARY OF PERTURBATION DATA 


VARIABLE NO. 
PERTURBED VALUE 
VARIABLE NO. 
PERTURBED VALUE 
VARIABLE NO. 
PERTURBED VALUE 
VARIABLE NO. 
PERTURBED VALUE 


1 MEAN VALUE 
0 . OOGOOD+OO 

2 MEAN VALUE 
0 . 000000+00 

3 MEAN VALUE 
0.000000*00 

4 MEAN VALUE 
0 . OOOOOD*OQ 


o.oooood*oo 

0 . 0O000D*O0 
0 . OOOOCD + OO 
0. 000000*00 


STANDARD DEVIATION 
STANDARD DEVIATION 
STANDARD DEVIATION 
STANDARD DEVIATION 


0. OOOOOD*00 
0 . 000000*00 
0 . 000000*00 
0.000000*00 


PERTURBATION ASSEMBLY. INCA. 0 PERT. 4 ITER. 0 CPTIKE- 


•** PERTURBATION 4 HAS NO EFFECT ON THIS ANALYSIS 

\ 


DATABASE UPDATE. INCR . 0 PERT. 4 ITER, 

1C . . . HEAT TRANSFER FOR A REGRIGERANT PIPE 
: 14- 5-199? 20:30 RESULTS PAGE: 1 


0 CPTIME" 


0.31 SEC 


0.32 SEC 


DATE 


TOTAL DISPLACEMENTS 


INCREMENT 0 TOTAL TRANSIENT TIKE 0.000000*0 


NODE 
COMP. 6 

COMP. 1 

1 

. OOOOOD+QO 

0 , OOOOOD+OO 

2 

. OOOOOD+OO 

0. OOOOOD+OO 

3 

.000000*00 

0. OOOOOD+OO 

4 

. 00C00D*00 

0. OOOOOD+OO 

5 

.000000*00 

0. OOOOOD+OO 

6 

. 00000D*D0 

0. OOOOOD+OO 

1C . . . HEAT 

TRANSFER FOR J 

> DATE: 14- 

5-1997 20:50 


COMP. 2 

COMP. 3 

0.70000D+02 

0. OOOOOD+OO 

0 . 35298D+02 

0.000000+00 

0. 352980+02 

0. OOOOOD+OO 

0 . 50614D+01 

0. OOOOOD+OO 

0 . S0614D+01 

0 . OOOOOD+OO 

0.5O009D+01 

0. OOOOOD+OO 


REGRIGERANT PIPE 
PAGE: 11 


COMP. 4 

COMP. S 


0. OOOOOD+OO 

0. OOOOOD+OO 

0 

0. OOOOOD+OO 

0. OOOOOD+OO 

0 

0. 000000+00 

0, OOOOOD+OO 

0 

0. OOOOOD+OO 

0. OOOOOD+OO 

0 

0. OOOOOD+OO 

0. OOOOOD+OO 

0 

0. OOOOOD+OO 

0. OOOOOD+OO 

0 


VERSION 6.01109 


END OF INCREMENT, INCR. 0 PERT. 0 ITER. 

1C . . . HEAT TRANSFER FOR A REGRIGERANT PIPE 
) DATE: 14- 5-1997 20:50 PAGE; 12 


0 CPTIME" 


0.34 SEC 

VERSION 6.0(109 


35 
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Appendix B 


FEM FILES FOR FLUID FLOW IN A DUCT 

(Refer to Chapter III, Fig. 18) 

NESSUS/FEM input file 


flowfaml.dat Sat Ju n 28 22:26:43 1997 


1 


•FEM 

C ... FLOW IN A CONSTANT DIAMETER PIPE {with A/L=0.1) 
C 
c 
c 

C . . . PARAMETER DATA 
C 

•BOUN 20 

•CONS 0 

•DISP 

•NODE 2 

*ELEH 1 

98 

•FORC 1 

•PERT 4 4 

•PRIN 

•BEAM 

•END 

€•••* MODEL DATA 


•COOA 

1 0.0000 

2 0.0000 

•ELEN 98 

1 1 2 

•BOUN 

2 1 0.0000 

2 3 0.0000 

1 l 0.0000 

1 2 1838.7 

1 3 0.0000 

1 4 0.0000 

1 5 0.0000 

1 6 0.0000 

•BEAM 1 

1 2 1.0000 

•ITERA 0 4 

40 0.0500 

•PROP 98 

1 2 1.0000 

•FORC 

2 2 -240. 2S 

•PRIN 

TOTA NCOE 


0.0000 

10.000 


1.0000 


14 89fi0 


0.0000 0.0000 

0.0000 0.0000 


0.3000 0.0000 


•END 


0.0000 

0.0000 


2 


. 5E- 4 


1.0000 

1 . 0000 


1 
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NESSUS/FEM output file 


£ 1 owfaml . out 


Sat Jtm 28 22:27:00 1997 


) DATE: 28- 6-1997 22:26 PAGE: 


9 


4 


•• PERTURBATION NO, 4 •• 


SUMMARY OF PERTURBATION DATA 


VARIABLE NO. 
PERTURBED VALUE 

1 KEAN VALUE 

0.000000*00 

0 .OOOOOD+OO 

STANDARD DEVIATION 

0. OOOOOD+OO 

VARIABLE NO. 
PERTURBED VALUE 

2 KEAN VALUE 

0 . OOOOOD+OO 

0. OOOOOD+OO 

STANDARD DEVIATION 

0. OOOOOD+OO 

VARIABLE NO. 
PERTURBED VALUE 

1 KEAN VALUE 

0.000000*00 

0. OOOOOD+OO 

STANDARD DEVIATION 

0 . OOOOOD+OO 

VARIABLE NO. 
PERTURBED VALUE 

4 KEAN VALUE 

0. 000000*00 

0. OOOOOD+OO 

STANDARD DEVIATION 

0 .000000+00 


PERTURBATION ASSEMBLY, INCR. 0 PERT. 4 ITER. 0 CPTIME. 
PERTURBATION 4 HAS NO EFFECT ON THIS ANALYSIS 


DATABASE UPDATE, INCR. 0 PERT. 4 ITER. 0 CPTIHE* 
1C ... FLOW IN A CONSTANT DIAMETER PIPE (with A/L«0.1) 

E: 28- 6-1997 22:26 RESULTS PAGE: 1 


0.27 SEC 


DAT 


00 


TOTAL DISPLACEMENTS 


INCREMENT 0 TOTAL TRANSIENT TIME Q.OOOOOD+ 


NODE COMP. 1 COMP 2 

COMP . 6 


1 0.000000*00 0 . 18397D+04 

0.000000*00 

2 0.000000*00 0.167740+04 

0. OOOOOD+OO 

IC ... rLOW IN A CONSTANT DIAMETER PIPE 
9> DATE: 28- 6-1997 22:26 PAGE: 10 


COMP. 3 

COMP . 4 

COMP. 5 

0 . OOOOOD+OO 

0. OOOOOD+OO 

0 . 000Q0D+00 

0. 000000+00 

0. OOOOOD+OO 

0. OOOOOD+OO 

eh a/l-0.1) 


VERSION 6.0(10 


END OF INCREMENT, INCR. 0 PERT. 0 ITER 0 CPTIME* 

C ... FLOW IN A CONSTANT DIAMETER PIPE (with A/L-0 1) 

> DATE: 2«- 6-1997 22:26 PAGE: U 


0.29 SEC 

VERSION 6.0(10 


45 


STOP DUE TO END OF INPUT FILE 
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Appendix C 


PFEM FILES FOR STRUCTURAL RELIABILITY 

(Chapter IV, Data is shown in Fig. 16 of Chapter III) 

NESSUS/PFEM input and output files for single Z-level 


•xp20pf «nw d&t Tu* Jun 24 18:13 i 00 1997 1 

•prat 

C ... PROBABILISTIC ANALYSIS FOR STRUCTURAL RELIABILITY OF EXP2 
•HVDBFINR 

•COND 0 

•DATATYPE 0 
•RESPTYFE 3 
•COMP 3 

•NODE 44 

•PERT 3 

1,2,3 

•RANVAA 3 

1,2,3 

•END 

•2FDKPINE 

•COMPUTATIONAL 1 3 

1,3,3 

•END 

•RVDEPINB 

•DEPINE 1 
PI 

1377.4 315.48 NORMAL 

PRESSURE 

30 5 6 1.0 

•DEFINE 2 

EMOD 

0 . 1700000E+08 0.17000008+07 NORMAL 

PROP 133 

1 58 0.0 1.0 0.0 0.0 0.0 

•DEFINE 3 

ALTA 

0 . 930Q000E-05 0 . 1900000E-05 NORMAL 

PROP 153 

1 38 0.0 0.0 0.0 1.0 0.0 
•PERT 1 
1 0.1 

•PERT 2 

2 0.1 

•PERT 3 
3 0.1 

•END 
•END 

c * ****** 

*F*M 

C . . .REFRIOERANT PIPE UNDER INSIDE PR E SURE 

C 

C 

c 

C . . . PARAMETER DATA 
C 

•CONS 0 
•DISP 

•NODE 58 
•EL EM 38 
153 

•BOUN 20 

•PRESSURE 

•TEMPERATURE 

♦LOUR 2 1/3 1 

•PERT 10 10 

PRESSURE 
PROPERTIY 
PROPERTY 
•PRIN 
•END 

C***' MODEL DATA 
•COOR 

1 -5,000 0.95 


2 

-4.500 

0.95 

3 

-4.000 

0.95 

4 

-3.500 

0.95 

5 

-3.000 

0.93 

4 

-2.700 

0.95 

7 

•2.400 

0.95 

8 

-2.100 

0.95 

9 

-1.800 

0.95 

10 

-1.500 

0.95 

11 

-1.200 

0.93 

12 

-0.900 

0.95 

13 

-0.600 

0.95 

14 

-0.300 

0.95 

15 

0.000 

0.93 

16 

0.300 

0.95 

17 

0.600 

0.95 

IS 

0.900 

0.95 

19 

1.200 

0.95 

20 

1.500 

0.95 

21 

1.800 

0.95 

22 

2.100 

0.95 

23 

2.400 

0.95 

24 

2.700 

0.95 

25 

3.000 

0.95 

26 

3.300 

0.95 

27 

4 . Oofe 

0.95 

28 

4.500 

0.93 

29 

5.000 

0.95 

30 

-5.000 

0.665 

31 

-4.500 

0.665 

32 

-4.000 

0.663 

33 

-3.500 

0.663 

34 

-3.000 

0.665 

35 

-3.700 

0.665 

36 

-2.400 

0.665 

37 

-2.100 

0.665 

38 

-1.800 

0.6(5 

39 

-1.500 

0.665 

40 

-1.200 

0.665 

41 

-0.900 

0.66S 

42 

-0.600 

0.665 

43 

-0.300 

0.663 

44 

0.000 

0.663 

43 

0.300 

0.665 

46 

0.600 

0.665 

47 

0.900 

0.663 

48 

1.200 

0.663 

49 

1.500 

0.665 

50 

1.800 

0.665 

31 

2.1O0 

0.665 

S2 

2.400 

0.665 

53 

2.700 

0.665 

54 

3.000 

0.665 

53 

3.500 

0.665 

56 

4.000 

0.665 

57 

4.500 

0.665 

58 

5.000 

0.(65 

ELEH 

153 


1 

1 30 

31 2 

2 

2 31 

32 3 

3 

3 32 

33 4 

4 

4 33 

34 5 

3 

5 34 

35 6 

6 

6 35 

36 7 

7 

7 36 

37 8 

8 

8 37 

38 9 

9 

9 38 

39 10 
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•xp2«pf am. dat 


Tua Jun 24 18:13:00 1997 


2 


10 

10 

39 

40 

11 

11 

11 

40 

41 

12 

13 

12 

41 

42 

13 

13 

13 

42 

43 

14 

14 

14 

43 

44 

15 

13 

13 

44 

43 

16 

1< 

16 

45 

46 

17 

17 

17 

46 

47 

18 

18 

18 

47 

48 

19 

19 

19 

48 

49 

30 

20 

20 

49 

50 

21 

31 

21 

50 

51 

22 

23 

22 

51 

52 

23 

23 

23 

52 

53 

24 

24 

24 

53 

54 

25 

25 

25 

54 

55 

26 

26 

26 

35 

56 

27 

27 

27 

56 

57 

28 

28 

28 

57 

58 

29 

'PRESSURE 





30 

58 

1577.4 



'TEMPERATURE 




1 29 

4.06 



30 58 

4.00 




•BOON 

IS 1 0.0 

44 1 0.0 

•ITERA 1 3 

40 0.0300 

•PROP 153 

1 58 0.1000B«-01 0 . 1700B+08 0.340 0.9500E-05 0.3188 6. 08*06 

•PR1N 

STRE NODE 

•■HD 

C *•*•• 

*PH 

HEAT EXCKAOER OP EXP2 WITH DATASETS 
•RVNUM 3 
•GFUNCTIOK 1 
•DATASETS 4 
•METHOD 1 
•PAIHTOPT 0 

•ANAL TYPE 1 

•END 

•ZLXWLS 1 
0.800 000*404 
•END 


DESIGN SENSITIVITIES 


TAYLOR SERIES EXPANSION OF THE PORK 

O « A0 ♦ A1*X1 ♦ A2»X2 ♦ ... ♦ AN*XN 
WHERE: 


axp2»pfaa.out Tua Jun 24 18:13:42 1997 3 

A0 IS THE CONSTANT TERM 

Al, A3, ... , AN ARE THE DESIGN SENSITIVITIES 
XI, X3, ... , XN ARE THE RANDOM VARIABLES 


TAYLOR SERIES COEFFICIENTS RANDOM EXPANSION NORN. DBS ION 

a co epps. value variable point <mpp» sjnsitivity 


0 0.3096481404 

1 0.3333341401 1 0.1S7740E+04 0.835343B4QO 
3 •0,1333311*03 3 C.170000E+08 -0.3S18S9E*00 
3 0.188750*407 3 0.930000E-05 0.276730E-02 


PERFORMING PROBABILISTIC ANALYSIS WITH FPI 


CDF RESULTS 

* U PROBABILITY ITER. NO. 

0 , 80000000E*04 0.56223765*401 0.99999999**00 0 
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NESSUS/PFEM input and output file for CDF 


exp2scdf.dat Tue Jun 24 19 j 19 *57 1997 1 

•ppm 

C ... PROBABILISTIC ANALYSIS FOR STRUCTURAL RELIABILITY 07 EXP2 
•KVD SPINS 

•COND 0 

•DATATYPE 0 
•RJESPTYPE 3 
•COMP 3 

•NODE 44 

•PERT 3 

1.2.3 

•AANVAR 3 

1.2.3 

•END 

•2PDE7INB 

•COMPUTATIONAL 1 3 

1,2.3 

•END 

•RVDEPINE 

•DEPINE 1 
PI 

0. 15774000E*04 0 . 3 154800B*03 NORMAL 

PRESSURE 

30 58 1.0 

•DEFINE 2 

EMOD 

0 . 1700Q00E+08 0. 1700000E + 07 NORMAL 

PROP 153 

1 58 0.0 1.0 0.0 0.0 0.0 

•DEFINE 3 

ALFA 

0, 9300000E-05 0 . 1900000E-05 NORMAL 

PROP 153 

1 31 0.0 0.0 0.0 1.0 0.0 

•PERT 1 

1 0.1 

•PERT 2 

2 0.1 

•PERT 3 

3 0.1 

•END 

•END 

c 

•FEM 

C . . .RTPRIGERANT PIPE UNDER INSIDE PRESURE 

C 

C 

c 

C , . . PARAMETER DATA 
C 

•CONS 0 
•DISP 

•NODE 38 
•ELEM 28 
133 

•BOUN 20 

•PRESSURE 

•TEMPERATURE 

•LOUB 2 1 • 3 1 

•PERT 10 10 

PRESSURB 
PROPKRTIY 
PROPERTY 
•PRIM 
•END 

C***» MODEL DATA 
•COOR 

1 '• -5.000 ~ 0.95 



2 

-4.500 

0.95 

3 

-4.000 

0.95 

4 

-3.500 

0.95 

5 

-3.000 

0.95 

6 

-2.700 

0.95 

7 

-2.400 

0.95 

8 

-2.100 

0.95 

9 

-1.800 

0.95 

10 

-1.500 

0.95 

11 

-1.200 

0.95 

12 

-0, 900 

0.95 

13 

-0.600 

0.95 

14 

-0.300 

0.95 

15 

0.000 

0.95 

16 

0.300 

0.95 

17 

0.600 

0.95 

11 

0.900 

0.95 

19 

1.200 

0.95 

20 

1.500 

0.95 

21 

1.800 

0.95 

22 

2.100 

0.95 

23 

2.400 

0.95 

24 

2.700 

0.95 

2S 

3.000 

0.95 

26 

3.500 

0.95 

27 

4 . 00(T 

0.95 

28 

4.500 

0.95 

29 

5.000 

0.9S 

30 

-5.000 

0.665 

31 

-4.500 

0.665 

32 

-4 . 000 

0.665 

33 

-3.500 

0.(65 

34 

-3.000 

0.665 

35 

-2.700 

0.6(5 

36 

-2.400 

0.665 

37 

-2.100 

0.665 

38 

-1.800 

0.665 

39 

-1.500 

0.665 

40 

-1.200 

0.665 

41 

-0.900 

0.665 

42 

-0.600 

0.665 

43 

-0.300 

0.665 

44 

0.000 

0.665 

45 

0.300 

0.665 

46 

0.600 

0.665 

47 

0.900 

0.665 

48 

1.200 

0. 665 

49 

1.500 

0.665 

50 

1.800 

0.665 

51 

2.100 

0.665 

52 

2.400 

0.665 

53 

2.700 

0.665 

54 

3.000 

0.665 

55 

3.500 

0.665 

56 

4.000 

0.665 

57 

4.500 

0.665 

58 

5.000 

0.665 

ELEM 

153 


1 

1 30 

31 2 

2 

2 31 

32 3 

3 

3 32 

33 4 

4 

4 33 

34 5 

5 

5 34 

35 6 

6 

6 35 

36 7 

7 

7 36 

37 8 

a 

e 37 

38 9 

9 

9 38 

39 10 
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2 


•xp2aedf.dat 


Tu« JMn 24 19 : 19 r 57 1997 


10 

10 

39 

40 

11 

11 

n 

40 

41 

12 

12 

12 

41 

42 

13 

13 

13 

42 

43 

14 

14 

14 

43 

44 

15 

13 

15 

44 

45 

16 

16 

16 

45 

46 

17 

17 

17 

46 

47 

18 

18 

18 

47 

48 

19 

19 

19 

48 

49 

20 

20 

20 

49 

30 

21 

21 

21 

50 

31 

22 

22 

22 

51 

52 

23 

23 

23 

52 

53 

24 

24 

24 

S3 

54 

25 

25 

25 

54 

55 

26 

26 

26 

35 

36 

27 

27 

27 

56 

57 

28 

28 

28 

57 

38 

29 


•PRESSURE 

30 58 0 . 1577400E+04 

•TEMPERATURE 

1 29 4.06 

30 58 4.00 

*BOUN 

15 1 0.0 

44 1 0.0 

•ITERA 1 3 

40 0.0500 

•PROP 153 

1 58 0.1000B+01 0 . 17OOE*08 0.340 0.9300E-05 

* PRXN 

STUB NODS 

•END 



*PPI 

HEAT EX CHAO EH OF EXP2 WITH DATASETS 
•RVNUM 3 
•OFUNCTION 1 
•DATASETS 4 
•METHOD 1 
•PHINTOPT 0 
•AMALTYPK 0 
•END 
•END 


_ m _* * f. r ° N s! NSITIV1TIbs 

TAYLOR SERIES EXPANSION OF THE FORM 
O - A0 ♦ A1*X1 ♦ A2*X2 ♦ ... ♦ wi*XN 
WHERE; 


•xp2«cdf . out Tue Jun 24 18 t 28 t 49 1997 

A0 IS THE CONSTANT TIM 

W| *3 an are the design sensitivities 

XI, X2 , ... , XN ARE THE RANDOM VARIABLES 


TAYLOR SERIES COEFFICIENTS 
A COBFFJ . VALUE 


RANDOM EXPANSION NORM. DESIGN 

VARIABLE POINT (KPP) SENSITIVITY 


0 0 . 20964SE+04 

1 0 . 2333348*01 1 
3 -0. 123321E-03 2 
3 0.168750E*07 3 


0 . 157740E*04 0.635343E*00 
0. 1700006*08 “0 . 3618B9E+00 
0.950000E-05 0. 276730E-02 


PERFORMING PROBABILISTIC ANALYSIS WITH FPI 


Z 

-0 . 206199418+03 
0. 55849989E+03 
0. 13238992E+04 
0. 2O8J29I5E+04 
0.285449788+04 
0. 43854964B+04 
0.5150S957E+04 
0 . 591 6293OE*04 
0.668169438*04 
0 . 74470935E+04 


CDF RESULTS 
U 

-0.510000008*01 
-0 . 41000000B+01 
-0.310000008*01 
-0.210000001*01 
-0.110000001*01 
0.900000008*00 
0 . 190000008*01 
0 . 29000000E+01 
0 , 390000006*01 
0. 49000000E+01 


PROBABILITY 
0 . 170122 3 lB-06 
0 . 20668716B-04 
0. 96767122E-03 
0.17864357E-01 
0 . 13 56661OE + O0 
0 . 81593991E*00 
0 . 97128351E+00 
0 , 996134 12E*00 
0 . 999951 88E*00 
0 . 99999952E* 00 


ITER. NO. 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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Appendix D 


PFEM FILES FOR THE HEAT TRANSFER RELIABILITY 

(Chapter IV, data is shown in Fig.16 in Chapter III) 

NESSUS/PFEM input and output files for single Z-level 


•xp2t*BJpf*mz.d*t Sun Jun 22 01:05i39 1997 1 

* PFEM 

q PROBABILISTIC ANALYSIS FOR TEMPERATURE FAILURE OP EXP 2 Z-LEVEL 

•MVDEPINE 

•COND 0 

♦DATATYPE 0 
•RESPTYPE 1 
♦COMP 2 

•NODS « 

•PERT 3 

1.2.3 

•RANVAR 3 

1.2.3 

•END 

•ZFDE7INS 

♦COMPUTATIONAL 1 3 

1.2,3 

♦«ND 

•RVDEPINE 

•DEFINE 1 

EMOD3 

0 . 39955476000E + 05 0 . 39955474000E+04 NORMAL 

PROP 98 

5 S 0.0 1.0 0.0 0.0 0.0 

•DEFINE 2 

EMOD2 

0 . 5887580008*02 0 . 588758000E+01 NORMAL 
PROP 98 

3 4 0.0 1.0 0.0 0.0 0.0 

•DSP INK 3 

ZMOD1 

0 . 785398E+02 0.785398K»01 NORMAL 

PROP 90 

i a o.o l.o o.o o.o o.o 

•PERT 1 

1 O.l 

•pert a 

2 0.1 

•PERT 3 
3 0.1 


♦END 

c * ••*••••••• 

*P*M 

C . . . HEAT TRANSFER FOR A REORIOERANT PIPE EXP2 

C 

C 

c 

c . . . PARAMETER DATA 
C 

•BOON 20 

•CONS 0 

•DISP 

♦NOD* « 

•DUPL 4 

•ELEN 1 

98 

•FORC 1 

♦PERT 10 10 - 

PROPERTIES 
PROPERTIES 
PROPERTIES 
•PRtN 
•BEAM 
•END 

C** •* MODEL DATA 
♦COOR 

1 0.0000 0.0000 0.0000 0.0000 0,0000 1.0000 
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2 

0.0000 

10.0000 

0.0000 

0.0000 

0.0000 

1.0000 

3 

0.0000 

10.0000 

0.0000 

0.0000 

0.0000 

1.0000 

4 

0.0000 

20.0000 

0.0000 

0.0000 

0.0000 

1.0000 

5 

0.0000 

20.0000 

0.0000 

0.0000 

0.0000 

1.0000 

6 

0.0000 

30.0000 

0.0000 

0.0000 

0 . 0000 

1 . 0000 


•ELEM 98 

1 X 2 

2 3 4 

3 5 6 

•DUPL 

3 2 

5 4 
’BOUN 

6 1 0.0000 

6 3 0.0000 

1 1 0.0000 

1 2 70.0000 

1 3 0.0000 

1 4 0.0000 

1 5 0.0000 

1 6 0.0000 

•BEAM 1 

1 6 1.0000 1.0000 

•ITEAA 1 3 

40 0.0500 

•PROP 98 \ 

1 2 1.0000 78.5398 0.3000 0.0000 2. St-4 

3 4 1.0000 58.8758 0.3000 0.0000 2. 52-4 

5 6 1.0000 39955.478 0.3000 0.0000 2.5E-4 

•FORC 

6 2 -241.68 

•PRIM 

TOTA NODE 

•END 

c * 

•FPI 

RBTRIGERANT PIPE RELIABILITY WITH Z- LEVEL (7 < 5«F1 
•RVNUW 3 
•GFUNCTION 1 
•DATASETS 4 
•METHOD 1 
•PR INTO PT 0 
•ANALTYPE 1 
•END 

*Z LEVELS l 

5.00000 

•END 


•xp2t*rapf •rnx .out Sun Jun 22 01t06il6 1997 3 

A0 IS THE CONSTANT TERM 

Al, A3, ... .AN ARE THE DESIGN SENSITIVITIES 
XI. X2 XN ARE THE RANDOM VARIABLES 


TAYLOR SERIES COEFFICIENTS RANDOM EXPANSION NORM. DESIGN 

A C08FFS. VALUE VARIABLE POINT {MPP> SENSITIVITY 


0 -0 . 7375068+02 

1 0. 151374E-OS 1 0 .3995558*05 0.841489E-03 

2 0 . 6971588*00 2 0.5887588*02 0.5710688*00 

3 0.3917658*00 3 0.7853988*02 0.4280908*00 


PERFORMING PROBABILISTIC ANALYSIS KITH FPI 


CDF RESULTS 

X U PROBABILITY ITER. NO. 

0.500000008*01 0 . 134027I3E*01 0.909921408*00 0 
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NESSUS/PFEM input and output files for CDF 


•xp 2 temp fern. dat 


Thu Apr 6 17tl4»4B 1995 


C PROBABILISTIC (COP) ANALYSIS FOR TEMPERATURE FAILURE OF EXP2 

•HVDEPIW1 

•COND 0 

•DATATYPE 0 
•RS3PTYPB 1 
•COHP 2 

•NODE 6 

•PERT 3 

1,2,3 

•RANVAR 3 
1,2,3 

•END 

•2FDEFINB 

•COMPUTATIONAL 1 3 

1.2,3 

•END 

•RVDBPINE 

•DEFINE 1 


EMOS3 

0.39955478000E*05 0 . 39955478000E+03 
PROP 98 

5 6 0.0 1.0 0.0 0.0 

•DEFINE 2 


NORMAL 

0.0 


EHOD2 

0.58675B000E*02 0 . S80758OOOE+O1 
PROP 98 

3 4 0.0 1,0 0.0 0 
•DEFINE 3 


NORMAL 

0 0.0 


EMOD1 

0 . 785398E*02 
PROP 98 

1 2 0.0 

• PERT 1 

1 0.1 

•PERT 2 

2 0.1 

•PERT 3 

3 0.1 


0 . 78S398E+01 NORMAL 

1.0 0.0 0.0 0.0 


•END 
•END 
C •• 


•PEN 

C . .. NEAT TRANSFER POR A RBGRIGERANT PIPE EXP2 

C 

c 

C 

C . . . PARAMETER DATA 
C 

•BOON 20 

•CONS 0 

•DISP 

•NODS « 

•DUPL 4 

*EL EM 3 

98 

•FORC 1 

•PERT 10 10 ' 

PROPERTIES 

PROPERTIES 

properties 

•PRIN 

•BEAM 

•END 

C*«* MODEL DATA 


•COOR 

1 


0.0000 0.0000 0.0000 


0.0000 0.0000 1.0000 


1 
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•ELEM 

1 

2 

2 

*DUPL 

3 

5 

*BOUN 

« 

6 
1 
1 


• ITERA 
40 

•PROP 

1 

3 

5 

•PORC 

6 

•PR IN 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 2 

3 i 

5 6 


0.0000 

0.0000 

0.0000 

■70.0000 

0.0000 

0.0000 

0.0000 

0.0000 


« 1.0000 

1 3 
0.0500 

9 8 

2 1.0000 

4 1.0000 

6 1.0000 


10.0000 

10.0000 

20.0000 

20.0000 

30.0000 


1.0000 


70.539 8 
58.1758 
39935.47. 


0. 0000 
0. 0000 
0.0000 
0.0000 
0.0000 


0.3000 

0.3000 

0.3000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0000 
0 0000 
0.0000 


0.0000 
0.0000 
0.0000 
0.0000 
0 . 0000 


2.5E-4 

2.51-4 

2.SE-4 


1 . 0000 
1.0000 
1.0000 
1.0000 
1.0000 


C **•••*•«»•«• 

•FPI 

HEAT EXCHAGER OP EXP2 WITH DATASETS 
•RVNUM 3 
*G FUNCTION l 
•DATASETS 4 
•METHOD 1 
•PRINTOPT 0 
’ANAL TYPE 0 
•END 
•END 


•xp2t*mpf am.out Mon May 19 22:40:3$ 1997 

A0 IS THE CONSTANT TERN 

*3. A2 AN ARE THE DESIGN SENSITIVITIES 

*2 XX ARE THE RANDOM VARIABLES 


TAYLOR SERIES COEFFICIENTS 
A COEFFS . VALUE 


RANDOH EXPANSION NORM. DESIGN 

VARIABLE POINT 1MPP) SENSITIVITY 


0 -0.7J7506E+02 

1 0 . 1S1374E-05 1 

2 0.5971581*00 2 

3 0.391765E*00 3 


0.399555E*0S 0.841489E-03 
0 . 588738E+02 0. 571068E+00 
O.705398E+O2 0 . 428090E*00 


PERFORMING PROBABILISTIC ANALYSIS WITH FPI 


Z 

-0.28039092£*02 
-0.229089411*02 
-0 . 177787908*02 
-0.12548639B+O2 
-0. 7518488 0E *01 
0 . 27418 138E*01 
0 . 787195 48E*01 
0.130021181*02 
0 . 181322678*02 
0.232524188*02 


CDP RESULTS 
U 

-0. 51000000E*01 
-0.41000000E*01 
-0. 310000001*01 
-0.21000000E*01 
-0. 11000000B*01 
0. 90000000E*00 
0.190000001*01 
0. 29000000E+01 
0 . 39000000E*01 
0.49000000E*01 


PROBABILITY 
0 . 17012231E-06 
0.206687HH-04 
0,967 67 122E-03 
0. 17864357E-01 
0. 135«6610B*00 
0. 81593991E*00 
0 . 97 12835 1S*00 
0 . 99813412E+00 
0 . 99995 1 B8E+00 
0 . 99999952E*00 


ITER. NO. 
0 
0 
0 
0 
O 
0 
0 
0 
0 
0 
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Appendix E 


PFEM FILES FOR RELIABILITY OF FLUID FLOW 


(Chapter IV, data is show in Fig.18 of Chapter III) 

NESSUS/PFEM input and output files for a single Z-level 


flowpf *mx . d&t Sat Jun 28 22t36t52 1997 

•PFEM 

C ... . PROBABILISTIC ANALYSIS FOR FLOW FAILURE 
•MVDEFINE 

•COND 0 

•DATATYPE 0 
•RESPTYPE 1 
•COMP 2 

•NODE 2 

•PERT 2 

1.2 

•RANVAR 2 

1.2 

'END 

•ZFDEFINE 

•COMPUTATIONAL 1 2 


•END 

•RVDEFINE 

•DEFINE 1 

EMOD 

14.096 1. 4896 NORMAL 

PROP 98 

1 2 0.0 1.0 

•DEFINE 2 

VELO 

-350 0 35.0 NORMAL 

FORC 

2 2 1.00 

•PERT 1 

1 0.1 
•PERT 2 
2 0.1 


0.0 


•END 

C 

•FEN 

C ... FLOW IN A CONSTANT DIAMETER PIPE (with A/L*0.11 

C 

c 

C V 

c . . . PARAMETER DATA 
C 

•BOUN 20 

•CONS 0 

•DISP 
•NODE 2 

•ELEM 1 

98 

•FORC 1 

•PERT 10 10 

VELOCITY 
FORCE 
•PR IN 
•BEAM 
•END 

C»»** MODEL DATA 

C00 J o 0000 0-0000 0.0000 0.0000 0.0000 1.0000 

3 0.0000 10.000 0.0000 0.0000 0.0000 1.0000 

•ELEM 98 

1 1 2 
•BOUN 

3 1 0 . 0000 

2 3 0.0000 

1 l 0.0000 

1 2 1838.7 

1 3 0.0000 
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1 4 0.0000 

1 5 0.0000 

1 6 0.0000 

•BEAM 1 

I 2 1.0000 1.0000 

•ITERA 1 2 

10 0.0500 

•PROP 98 

1 2 1.0000 14.8960 0.3000 0.0000 2 5E-4 

•FORC 

2 2 -350.0 

*PRIN 

TOTA NODE 

•END 

c . •••••• 

‘PPI 

FLOW FAILURE PROBABILITY WITH Z-LEVEL I P< 1677 . 4d*1 ) 
•RVNUM 2 
•OFUNCTION 1 

•DATASETS 3 

•METHOD 1 

•PRINTOPT 0 
•ANALYTYPE I 

•END 

*Z LEVELS 1 

1677.4 

•END 


DESIGN SENSITIVITIES 


TAYLOR SERIES EXPANSION OF THE FORM 

G * A0 ♦ Al'Xl ♦ A2 *X2 * ... ♦ AN'XN 
WHERE: 

A0 IS THE CONSTANT TERM 

Al. A2 AN ARE THE DESIGN SENSITIVITIES 

XI, X2 XN ARE THE RANDOM VARIABLES 


TAYLOR SERIES COEFFICIENTS RANDOM EXPANSION NORM. DESIGN 

A COEFF5. VALUE VARIABLE POINT (MPP) SENSITIVITY 


0 0.160374E*04 


flowpfamx .out flat Jun 28 22i37tl3 1997 3 

1 0 . 15773SE»02 I 0.148960E*02 0.5000006*00 

2 0.671321E*00 2 -0 . 150000E*03 -0 . 500000E*00 


PERFORMING PROBABILISTIC ANALYSIS WITH PPI 


CDF RESULTS 

Z U PROBABILITY ITER. NO. 

0.16774000E*04 O . 22 174074E*01 0 . 9B670241E*00 0 
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NESSUS/PFEM input and outpout files for CDF 


flowcdf.dat Sat Jun 28 22:45*21 1997 1 

*PFEM 

C PROBABILISTIC ANALYSIS FOR FLOW FAILURE 

•KVDEFINE 

•COND 0 

•DATATYPE 0 
•AESPTYPE 1 
•COMP 2 

•NODE 2 

•PERT 2 

1,2 

•RANVAR 2 
1.2 
•END 
•ZFDEFINE 

•COMPUTATIONAL 1 2 

1,2 

•END 

•RVDEFINE 

•DEFINE 1 

EXOD 

14.896 1.4896 NORMAL 

PROP 98 

1 2 0.0 1.0 0.0 0.0 0.0 

•DEFINE 2 

VELO 

-350.0 35.0 NORMAL 

FORC 

2 2 1.00 
•PERT 1 
1 0.1 
•PERT 2 
2 0.1 
* END 
•END 

C 

*PEM 

C ... FLOW IN A CONSTANT DIAMETER PIPE (with A/L=0.1t 
C 
C 
C 

C . . . PARAMETER DATA 
C 

•BOON 30 
•CONS 0 
•DISP 

•NODE 2 
•ELEN 1 
98 

•FORC 1 
•PERT 10 10 

VELOCITY 
FORCE 
•PRIN 
•BEAM 
•END 

€•♦•• HOOEL DATA 
•COOR 

1 0.0000 0.0000 0.0000 0.0000 O.OOOG 1.0000 

2 0.0000 10.000 0,0000 0.0000 0,0000 1.0000 

•ELZM 98 

1 1 2 
•BOUN 

2 1 0.0000 

2 3 0.0000 

1 1 0.0000 

1 2 1830.7 

l 3 0.0000 
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0.0000 

0.0000 

0.0000 


’BEAM 1 

1 2 1.0000 1 0000 
•ITEAA 1 2 

40 0.0500 

•PROP 98 

1 2 1.0000 14.896 0.3000 0.0000 2 5E-4 

*PORC 

2 2 -350.0 

•PRIN 

TOTA NODE 

'END 

C* •*•■*• 

*FPI 

FLOW FAILURE PROBABILITY WITH A SERIES 2 -LEVELS 
•RVNUM 2 
•GFUNCTION I 
•DATASETS 3 

•METHOD 1 

•PRINTOPT 0 
■AKALYTYPE 0 
•END 
•END 


DESIGN SENSITIVITIES 

TAYLOR SERIES EXPANSION OF THE FORM 

G - A0 * A1*X1 ♦ A2 *X2 ♦ ... * AN a XN 
WHERE ; 

A0 IS THE CONSTANT TERN 

Al, A2 , ... , AN ARE THE DESIGN SENSITIVITIES 
XI. X2 XN ARE THE RANDOM VARIABLES 


TAYLOR SERIES COEFFICIQ4TS 
A C0IFFS. VALUE 


0 0. L60374C.04 


RANDOM EXPANSION NORM. DESIGN 

VARIABLE POINT (MPP> SENSITIVITY 


flowcdf .out S* t Jun 

1 0.157735E+02 1 

2 0.«71321E*00 2 


28 22 i 46:36 1997 

0. 14B960H+02 0 . 5000COE»00 

-0.350000E+Q3 -0. SOOOOOE+OO 


PERFORMING PROBABILISTIC ANALYSIS WITH FPI 


Z 

0. 143431<«B*04 
0.146753652*04 
0. 15007S64S+04 
0. 15339763E + 04 
0. 15671 961E*04 
0. 1633635JE*04 
0. 1666B5S8E*04 
0. 17000756E*04 
0. 17332955E*04 
0. J7«S1S4E*04 


CDP RESULTS 
U 

-0 . 5X0000001*01 
-0.4 1000000E*01 
-0.3 lOOOOOOE+Ol 
-0.2 10000001+01 
-0 . 11000000E*01 
0. 90000000E*00 
0 . 19000000E*01 
0 . 2 >OOOOOOE*01 
0 . 39000000E+01 
0 . 49000000E+0 1 


PROBABILITY 
0. 17012231E-06 
0 . 206687161-04 
0.96767I22E-OJ 
0 . 17864 J57E-01 
0. 13566610E+00 
0 . 8 15 93991E*00 
0.971283518*00 
0. 998134128*00 
0. 99995108E*OO 
0. 99999952E+00 


ITER. NO. 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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Appendix F 


GFUN.DAT AND GFUN.MOV FILES OF IMPOSITIONS 

(Chapter IV, data is show in Fig. 16 of Chapter III) 

NESSUS/PFEM input gfun.dat and output gfun.mov files 



10 

- 1 .500 

0.95 

11 

- 1.200 

0. 95 

12 

- 0.900 

0.95 

13 

- 0.600 

0.95 

14 

- 0.300 

0.95 

15 

0 . 000 

0.95 

16 

0.300 

0.95 

17 

0.600 

0.95 

18 

0.900 

0.95 

19 

1.200 

0.95 

20 

1.500 

0.95 

21 

1.800 

0 . 95 

22 

2.100 

0.95 

23 

2.400 

0.95 

24 

2.700 

0.95 

25 

3 .000 

0. 95 

26 

3.500 

0.95 

27 

4.000 

0.95 

28 

4.500 

0.95 

29 

5.000 

0.95 

30 

- 5.000 

0.665 

31 

- 4.500 

0.665 

32 

- 4 . 000 

0.665 

33 

- 3 . 500 

0 . 665 

34 

- 3.000 

0.665 

35 

- 2.700 

0.665 

36 

- 2.400 

0.665 

37 

- 2.100 

0.665 

38 

- 1 .800 

0.665 

39 

-i.soo 

0. 665 

40 

- 1.200 

0.665 

41 

- 0.900 

0.665 

42 

- 0.600 

0.665 

43 

- 0.300 

0.665 

44 

0.000 

0.665 

45 

0.300 

0.665 

46 

0.600 

0.665 

47 

0.900 

0.665 

4 B 

1 .200 

0.665 

49 

1.500 

0.665 

50 

1.800 

0.665 

51 

2.100 

0.665 

52 

2.400 

0.665 

53 

2.700 

0.665 

54 

3.000 

0.665 

55 

3.500 

0.665 

56 

4.000 

0.665 

57 

4.500 

0.665 

58 

5 .000 

0.665 

'EL EM 

1 S 3 


1 

1 30 

31 2 

3 

2 31 

32 3 

3 

3 32 

33 4 

4 

4 33 

34 5 

5 

5 34 

35 6 

6 

6 35 

36 7 

7 

7 36 

37 8 

8 

8 37 

38 9 

9 

9 38 

39 10 

10 

10 39 

40 11 

11 

11 40 

41 12 

12 

12 41 

42 13 

13 

13 42 

43 14 

14 

14 43 

44 15 

15 

15 44 

45 16 

16 

16 45 

4 6 17 

17 

17 46 

47 18 
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2 


0furk.dat Turn 


18 


18 

47 

48 

19 


19 

48 

49 

20 


20 

49 

50 

21 


21 

50 

31 

22 


22 

51 

32 

23 


23 

52 

53 

24 


24 

53 

54 

25 


25 

54 

55 

26 


26 

55 

56 

27 


27 

56 

57 

28 


28 

37 

58 

•PRESSURE 




30 

58 


3300.0 


•TEMPERATURE 




1 

29 


4.06 


30 

50 


4.0 


•BOUN 





15 

1 

0 

.0 


44 

1 

0 

.0 


•ITERA 

1 

2 



40 

0.0500 



•PROP 

153 




1 

58 

0 

. 100 08*0 1 


•PRIM 

STUB NODE 


Jul 1 15 i 39 i 34 1997 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


0 . 17008+08 0.340 0.95008-05 


'END 
C *•< 
•FPI 


HEAT EXCHAOER OF EXP2 KITH DATASETS 
•RVNUM 2 
* a FUNCTION 1 
•DATASETS 3 
•METHOD 1 
•PRINTOUT 0 
•ANAL TYPE 1 
•END 

•ZLBV8LS 1 
0 . 8000008*04 
•END 


0.3X68 


6 . 0E + 06 


gfun.&ov 

CONDITION « 
I NCR 0 


Tu« Jul 1 15 1 28 1 37 1997 

NODE COMPONENT LAYER 

“ j V 


1 


0.6549378*04 0.6549378*04 0.2125248*03 


NUMBER or PROBABILITY LEVELS 

ZO. PROS. - 0.8000008*04 1.000000000 

JMOD 0.53976318*07 0.999887 -6.824923 

***** 0.9695273B-05 -0.015037 0.102775 


CDF RESULTS 
z U 

0.800000008*04 0.682569568*01 


PROBABILITY 

0.100000008*01 


ITER. NO. 
0 
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gUun.mov 


Sun Jun 29 03i16i45 1997 


1 


CONDITION 
I NCR 


COMPONENT 

3 


LAYER 

I 


Z MEDIAN, KEAN, STD. DEV. 

0 . 678271E+04 0.678271E+04 0 . 21 2754E+03 

NUMBER OP PROBABILITY LEVELS 1 

LVL, 20, u, PROB. - 1 0 - 800000E+04 5.72159 0.999999995 

DESIGN POINT COORD., SENSITIVITY FACTORS, STD. DEV. FROM MEAN VALUES 
ZD, PROB. • 0 . 800000E*Q4 0.999999995 

EMOD 0 .7274407E+07 0.999886 -5.720937 

ALFA 0.9664020B-05 -0.015088 0.086326 


CDF RESULTS 
Z U 

0 .800000002*04 0 . 57215881E+01 


PROBABILITY 

0.999999998*00 


gfun.mov 


Bun Jun 29 03*llt36 1997 


CONDITION « NODE 


COMPONENT LAYER 


1NCR 0 44 3 


1 


Z MEDIAN, MEAN, STD. DEV. 

0 . 701604E*04 0.701604E+Q4 0.212994E+03 

NUMBER OF PROBABILITY LEVELS 1 

LVL. Z0, u, PROB. • 1 0 . 800000E+Q4 4.61965 

DESIGN POINT COORD., SENSITIVITY FACTORS, STD. DEV. 
Z0, PROB. = 0.8000008*04 0.999990076 

BMOD 0.9147492E*07 0.999886 -4.619122 

ALFA 0 . 9632282E-05 -0.015071 0.069622 


0.999998076 
FROH MEAN VALUES 


Z 

0 . 80000000E*04 


CDF RESULTS 
U 

0.461964718*01 


PROBABILITY 

0.999998088*00 


gfun.mov Bun Jun 29 03il3sSO 1997 1 


CONDITION > NODS COMPONENT LAYER 


INCR 0 44 3 1 


2 MEDIAN. MEAN. STD. DEV. 

0 . 74827 1E*04 0.748271E*04 0.213454E*03 

NUMBER OF PROBABILITY LEVELS 1 

LVL. Z0. u, PROB. * 1 0.8000008+04 2.42343 0.992312628 

DESIGN POINT COORD., SENSITIVITY FACTORS, STD. DEV. FROM MEAN VALUES 

20, PROB. - 0 . 800000E+04 0.992312628 

EMOD 0 . 1288064E+09 0.999887 -2.42315S 

ALFA 0 . 9569244E-05 -0.015038 0.036444 


CDF RESULTS 

Z U PROBABILITY ITER. NO. 

0. BOOOOOOOE+04 0 . 24234292E+01 0.992312638+00 0 
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gf un.mov 


Turn Jul l 15 i 39 i 56 1997 


1 


CONDITION 

INCR 


COMPONENT 

3 


2 MEDIAN. MEAN, STD. DEV. 

0.7716048+04 0.7716048+04 

NUMBER OP PROBABILITY LEVELS 
LVL, ZO, u, PROB. a 1 0. 

DESIGN POINT COORD., SENSITIVITY FACTORS. STD 
20. PROB. - 0 . 000000004 0.908043163 

0400 0.14741298+08 0.999887 -1.328651 

ALFA 0.95379258-05 -0.015021 0.019961 


0.2135948+03 
1 

8000008+04 1.32880 

DEV. 


0.908043163 
FROM MEAN VALUES 


CDF RESULTS 
Z U 

0.800000008*04 0.132880121*01 


PROBABILITY 

0.908043168*00 


gf un . mov 

CONDITION I 
INCR 0 


Turn JUl 1 15:36:47 1997 

node component layer 

44 3 r 


1 


0 . 2139248+03 


0.7949388+04 0.794938E+04 

NUMBER OF probability LEVELS 

U ' PR0B- * 1 0 8000006+04 0.23664 

DESIGN POINT COORD., SENSITIVITY FACTORS, STD. DEV. 
20. PROS. ■ 0 . 800000B+04 0.593531990 

»OD 0.16597768+01 0.999880 -0.236614 

ALFA 0 . 9506726E-05 -0.014959 0.003540 


0.593531990 
FROM MEAN VALUES 


2 

0.800000008+04 


CDr RESULTS 
U 

0.236640198+00 


PROBABILITY 

0.593S31998+00 


ITRR. NO. 
0 


gfun.mov Tu« Jul 1 16:17:30 1997 1 

CONDITION 8 NODE COMPONENT LAYER 


INCR 0 44 3 


Z MEDIAN, MEAN, STD. DEV. 

0. 7996048+04 0.7996048+04 0.2139648+03 

NUMBER OF PROBABILITY LEVELS 1 

LVL, Z0, u, PROB. ■ 1 0 . 800000E+04 0.01849 0.507377567 

D83KW POINT COORD., SENSITIVITY FACTORS, STD. DEV. FROM MEAN VALUES 
Z0, PROB. - 0.8000008+04 0.507377567 

*M° D 0.16968568+08 0.999887 -0.018492 

ALFA 0.95005278-05 -0.015003 0.000277 


CDF RESULTS 

Z U PROBABILITY ITRR. NO 

0.80000000E+04 0 . 18493758H-0I 0.507377578+00 0 
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Appendix G 


DATA AND OUTPUT FILES FOR CALCULATION OF UNION 

PROBABILITY 


(Chapter IV) 


temp.dat 


Tue Jul 1 16:23: 17 1997 


1 , 2 , 1 

1, 17000000.0,1700000.0 

2, 0.0000095,0.0000019 

1, 6.825696 

2, 5.721588 

3, 4.619647 

4, 2.423429 

5, 1.328801 

6, 0.236641 

7 , 0.018497 

1, 6.825696 0.999887 -0.015057 

2, 5.721588 0.999886 -0.015088 

3 4 619647 0.999886 -0.015071 

4, 2.423429 0.999887 -0.015038 

5^ 1 328801 0.999887 -0.015021 

6, 0.236641 0.999888 -0.014959 

7, 0.018497 0.999887 -0.015003 


1 


********** BOUNDS ***** 
Lower bound = 0.4926211 
Upper bound = 0.4926211 
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